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ABSTRACT 

Exi  sti  ng  ai  rworthi  ness  standards  requi  re  that  al  I  aspects  of  fati  gue  be  consi  dered  i  n  the  desi  gn, 
development  and  maintenance  of  aircraft  structures.  To  minimise  costs  and  improve  combat 
readiness  requires  accurate  and  efficient  fatigue  life  evaluation.  This  report  details  the 
introduction  of  a  strain-life  algorithm  known  asFAMSH  into  the  software  tool  suite  known  as 
CGAP  devel  oped  and  mai  ntai  ned  at  the  DSTO .  A I  ong  w  i  th  i  mprovements  to  the  ori  gi  nal  FA  M  SFI 
code  this  most  recent  development  introduces  a  materials  database  and  dynamic  memory 
allocation.  With  other  additional  improvementsthelatestreleaseoftheCGAPenvironmentseeks 
to  improvethe  ease  with  which  theengineerisableto  accurately  and  reliablyconductfatiguelife 
analysis.  Contained  within  this  report  is  an  explanation  of  the  strain-life  theory  in  conjunction 
with  adetailed  description  of  its  implementation  intheCGAP  GUI  environment.  Examples  are 
provided  to  guidethe  user  through  its  operation. 
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Executive  Summary 

Theassessment  of  fatigue  damage  in  structures  is  an  essentiai  and  an  integrai  part  of  the 
management  of  air  vehicies  within  the  Austraiian  Defence  Force  (ADF).  The  ADF 
Austraiian  Air  Pubiication  7001.053  states  "...  to  assure  structurai  airworthiness, 
management  of  Aircraft  Structurai  i  ntegrity  (ASi )  is  essentiai  from  the  acquisition  stage, 
when  the  specification  is  deveioped,  through  to  fieet  retirement."  As  such,  the  aircraft 
certification  process  requires  that  fatigue  damage  be  assessed  and  criticai  areas, 
susceptibieto  fatigue  damage,  identified.  This  highiights  the  importance  of  methods  to 
evaiuatefatiguedamagetotheoveraii  safe  and  economicai  managementof  aircraft.  FI  ence 
theneedfor  robust  user  friendly  softwarewhich  theengineercan  use  to  assess  thefatigue 
life  of  critical  structural  elements. 

Improvements  in  these  methods  and  tool  scan  be  achieved  in  several  ways,  for  example, 
by  providing  consistent  material  data,  the  use  of  built-in  error  checking,  and  the  use  of  a 
fami  I  i  ar  graphi  cal  user  i  nterface.  This  docu ment  d i scusses the  i  ncorporati  on  of  a  strai  n-l  ife 
methodol  ogy  i  nto  an  exi  sti  ng  software  tool  CGA  P.  A  detai  I  ed  descri  pti  on  of  thestrai  n-l  ife 
theory  and  how  it  has  been  implemented  intotheCGAP  environment  is  provided.  The 
document  also  provides  a  thorough  explanation  of  the  operation  of  the  software  tool. 
Fi  nal  ly  a  step  by  step  probi  em  is  worked  through  for  the  users  benefit. 

A  ecu  rate  fati  gue  I  i  fe  assessment  and  I  ife  extensi  on  of  metal  I  i  c  ai  rcraft  structu  res  d  i  recti  y 
benefits  defence.  These  benefits  through  improvements  in  prediction  reliability  can  be 
realised  in  terms  of  reduced  operational  and  maintenance  costs,  improved  performance 
and  combat  readiness. 
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1.  Introduction 

The  assessment  and  assu  ranee  of  fati  gu  e  d  amage  pred  i  cti  ons  i  n  stru  ctu  res  i  s  an  essenti  al  part 
of  the  management  of  ai  r  vehi  cl  es  w  i  thi  n  the  A  ustral  i  an  Defence  Force.  The  A  A  P  7001.053  [  1] 
states  that  "to  assure  structural  airworthiness,  management  of  Aircraft  Structural  Integrity 
(ASI)  Isessentlal  fromtheacquisitlon  stage,  when  the  specification  Isdeveloped,  through  to 
fleet  retirement."  Thethrough-llfe  support  of  aircraft  necessitates  the  assessment  of  fatigue 
damage  caused  by  service  loads  and  Its  effects  on  structural  airworthiness.  As  such  methods 
of  fatigueevaluatlon  are  Important  In  the  overall  safe  and  economical  management  of  aircraft. 
In  practi  ce  the  analytical  assessment  of  fatigue  damage,  substantiated  by  coupon,  component 
and  full-scaletests,  form  the  basisof  technical  data  for  the  design  of  airworthiness  structural 
IntegrItymanagementplans.GIventhat  accu  rate  f  atl  g  u  e  1 1  f  e  assessment  I  s  cr  I  tl  cal  to  th  e  safe 
and  economical  through-llfesupport  of  aircraft.  It  Is  Important  that  we  develop  algorithms 
and  toolsthatareintultiveand  user-fri  endly  to  support  thestructural  engineer  In  the  process 
of  quality  assurance. 

The  total  fatigue  life  of  a  structure  may  be  assessed  by  using  oneof  four  approaches:  (1)  by 
assuming  the  total  fatigue  life  of  a  structure  Is  governed  entirely  by  crack  Initiation;  (2)  by 
assuming  that  the  total  fatigue  life  Is  governed  entirely  by  crack  growth;  (3)  by  treating  the 
total  fatigue  life  of  a  structure  as  a  serial  combination  of  crack  Initiation  and  crack  growth;  and 
finally  (4)  by  a  parallel  application!  of  both  crack  Initiation  and  crack  growth  analyses.  The 
firstapproach  wastypicallyused  Intheearlyyearsoffatigueassessmentand  Isstlll  thedesign 
practice  of  some  Industries  where  no  regular  Inspections  of  thestructuresarecarried  out.  The 
second  approach,  which  Is  an  emerging  method,  relies  on  the  accurate  prediction  of  the 
growth  of  very  smal  I  cracks  that  takes  pl  ace  I  n  the  early  stages  of  fatl  gue  damage.  With  the 
gradual  Improvement  In  short  crack  growth  theory,  this  approach  may  gain  more 
prominence.  The  serial  combination  approach,  which  has  been  adopted  for  the  P-3C  fleet 
assessment  [2],  treatstheearly  stages  of  fatigue  damage  as  a  process  of  crack  Initiation  and 
then  treats  thesubsequent  fatl  gue  damage  In  terms  of  crack  growth2.Thefourth  approach  Is 
being  used  on  some  later  generation  combat  aircraft,  particularly  where  low  structural 
maintenance  must  be  assured. 

Different  fatigue  crack  Initiation  models  are  routinely  used  In  the  safe-life  design  and  the 
development  of  aircraft  structural  Integrity  management  plans.  One  of  the  most  commonly 
used  fatigue  crack  Initiation  models  Is  the  so-called  strain-life  method.  The  theory  that 
underpins  the  strain-life  method  forms  the  basis  of  the  fatigue-llfe  tool  FA  MS  [3,  4]  and 
FA  M  SFI .  FA  M  SFI  Isa  mod  I  f  I  ed  versi  on  of  FA  M  S  w  I  th  ad  d  I  tl  onal  expand  ed  f  u  ncti  onal  I  ty  w  I  th 
fewer  1 1  mitatl  ons  I  mposed .  FI  owever,  both  FA  M  S  and  FA  M  SFI  I  ack  an  I  ntultl  ve  user  I  nterface 
and  amethod  by  which  to  manage  material  properties.  RecentlytheFAMSFI  sourcecodehas 
been  Incorporated  as  a  module  Into  the  Windows  GUI  known  asCGAP.  CGAP  provides  a 
common  user  friendly  utility  that  can  be  used  to  control  the  material  property  definitions 


!  The  particulars  of  the  combination  account  for  the  limitations  of  each  model  assessed  against  part 
Inspectablllty;  damage  susceptibility;  N  Dl  capability;  criticality  of  the  structure;  new  aircraft  materials; 
and  cost  to  repair/  replace 

2  The  crack  I  nitl  atl  on  and  growth  stages  are  del  I  neated  by  a  crack  size  relatl  ng  to  the  capabi  I  Ity  of  non¬ 
destructive  Inspection  technique  and  the  limitations  of  the  crack  growth  models. 
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which  are  important  in  achieving  consistency  in  input  data  that  underpins  the  advice 
provided  to  the  ADF.  The  introduction  of  error  checking,  consistent  materiai  database  and 
intuitive  interface  wiii  uitimateiy  minimise  mistakes  and  improve  the  accuracy  and 
assuredness  of  advice. 

The  ai  ms  of  th  i  s  d  ocu  ment  are  to 

•  summarise  the  saiient  points  of  the  strain-iife  theory; 

•  expiain  how  the  strain-iife  approach  has  been  impiemented  in  CGAP; 

•  provide  guidance  for  the  new  moduie; 

•  provideexampies  demonstrating  the  use  of  the  new  moduie. 

Accurate  fatigue  iife  assessment  and  iife  extension  of  metaiiic  aircraft  structures  directiy 
benefits  defence.  These  benefits  can  be  reaiised  in  terms  of  reduced  operationai  and 
maintenance  costs,  improved  performance  and  combat  readiness,  improvements  in  these 
methods  and  toois  can  be  achieved  in  severai  ways,  for  exampie,  by  providing  consistent 
materi ai  data,  bui  it-i  n  error  checki  ng,  and  a  fami  i  i ar  graphi  cai  user  i  nterface  (GU  i ).  A  i  i  these 
and  others  w  i  i  i  i  ead  to  i  mproved  robustness  i  n  pred i  cti  ons  w h i  ch  u  i  ti  matei  y  contri  butes  to  the 
more  efficient  and  reiiabie  management  of  A  Si  fatigue  advice. 

The  advice,  guidance  and  deveiopments  contained  in  this  report  contribute  directiy  and/  or 
indirectiy  to  thefoiiowing  important  high  ievei  issues: 

1.  Quaiity  assurance  is  an  integrai  part  of  ASi  management 

2.  A  wide  range  of  design  phiiosophies,  ASi  management  phiiosophies,fatiguemodeis, 
materiai  types  and  anaiysis  options  must  be  ad  dressed  tosupportADF  aircraft,  both 
in-being  and  in-acquisition. 

3.  Austraiiandefenceaircraftaresustained  by  muitipiepartners,  nameiytheindividuai 
OEMsand  theAustraiian  Defence  Support  Network  (DSN)^ 

4.  A  common, fiexibieandtransparentfatiguepredictiontooi  setisrequiredtoassureaii 
of  the  above  can  be  supported. 


2.  C G  A  P;  C  rack  G  rowth  A  nalysis  Program 

CGAPisaDSTO  devei  oped  softwaretooi  that  contai  ns  a  set  of  ai  gorithms  commoniy  used  to 
assess  and  anaiyse  fatigue  damage,  itisenvisaged  thatCGAP  wiii  graduaiiy  becomeaone- 
stop  tooi  set  with  which  the  engineers  can  perform  most  of  their  required  fatigue  anaiyses. 
CGAP  is  a  Microsoft  Windows- based  appiication  withaGUi  and  an  integrated  database  for 
the  management  of  materiai  properties,  geometry  sets  and  ioad  cases.  TheGUi  assists  the  user 
in  entering  input  data,  inspecting  spectra  and  checking  for  data  consistency,  itai  so  provides  a 
simpie  piotting  capabiiity  for  crack  growth  outputs.  CGAP  contains  a  native  crack  growth 
anai  ysi  s  mod  ui  e  based  on  the  concept  of  pi  asti  ci  ty-i  nd  uced  crack  ci  osu  re.  i  n  add  i  ti  on  to  ai  i  the 
functi  onai  i  ties  of  FASTRA  N  3.8[5],  the  native  moduie  aiso  hasthecapabi  i  ity  for  crack  growth 
anaiysis  invoiving  notch  piasticity,  and  probabiiistic  crack  growth  anaiysis  based  on  the 


3TheDSN  comprises  ASI -DGTA,  DSTO,  SPOsand  numerous  Australian  industry  partners 
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Monte  Carlo  method.  For  probabilistic  analysis,  the  initial  crack  size,  crack  growth  rate 
parameters  and  the  peak  spectrum  stress  can  be  randomised  using  three  distribution 
functions.  Importantly,  CGAP  interfaces  seamlessly  with  FASTRAN3.8[5]  and  theDKEFF[5] 
code.  A  recent  update  to  the  program  includes  the  incorporation  of  the  FAMSFI  code 
described  in  this  report.  The  inclusion  of  FAMSFI  to  the  CGAP  environment  expands  the 
functi onal  ity  of  CGA P  to  i ncl ude  a  strai  n-l  ife  approach  for  crack  i niti ati on  analysis. 


3.  FAMSH;  Strain  Life  Algorithm 

FAMSFI  is  a  software  tool  based  on  FAMS[3,  4]  and  developed  attheDSTO.  It  is  used  to 
perform  fatigue  life  calculations  based  on  a  strai  n-l  ife  methodology.  FAMSFI  was  originally 
developed  to  support  the  technical  interpretation  and  development  of  the  structural 
management  plan  (SM  P)  for  theP-3C.  It  is  currently  being  used  in  the  certification  process  of 
the  C-130J.  FAMSFI  is  based  on  FAMS[3, 4],  with  several  added  equivalent  strain  equations, 
new  outputs  for  damage  analysis  and  support  for  the  DBI/  SST  spectrum  format  [6].  The 
cu  rrent  rel  ease  of  FA  M  SFI  provi  d  es  f  u  1 1  y  dynami  c  memory  management  to  al  I  ow  the  program 
to  be  used  on  spectra  of  any  size,  only  limited  by  the  size  of  the  physical  memory  of  the 
hardware,  without  recompilation.  Furthermore,  thesameFAM  SFI  sourcecodeisnow  used  to 
generate  the  executable  file  on  both  Windows  and  Linux  platforms,  thus  ensuring  the 
consistency  and  currency  of  the  results  obtained  from  the  code. 


4.  The  Strain-Life  Approach 

Trad  i  ti  onal  I  y,  most  theori  es  on  fati  gue  rel  ate  the  d  amage  caused  by  cycl  i  c  I  oad  i  ng  to  the  I  ocal 
stress  or  strai  n  ranges,  result!  ng  i  n  two  si  mi  I  ar  but  different  approaches:  stress-1  ife  approach 
and  strai  n-l  ife  approach.  The  stress-1  ife  approach  is  commonly  used  forlow  loads  and  high 
cycle  fati  gue  when  local  plasticity  is  negligible,  while  the  strai  n-l  ife  approach  is  commonly 
used  in  situations  where  local  plasticity  may  occur.  For  the  stress-life  approach,  the  main 
material  data  required  is  the  stress-life  curve  for  the  stress  concentrator  concerned.  Thus 
different  stress-1  ife  data  need  to  be  generated  for  different  notch  features.  The  I  ocal  stresses 
and  strains  are  calculated  using  the  stress  concentration  factor  of  the  notch  feature.  For  the 
strai  n-l  ife  approach,  it  is  assumed  that  as  far  as  fati  gue  damage  i  s  concerned  the  materi  al  at 
the  notch  root  behaves  I  ike  as  a  smooth  (un-notched)  specimen.  Thus,  only  oneset  of  strain- 
life  data  obtained  from  a  smooth  specimen  is  required  for  the  basic  material  data.  Flowever, 
validation  against  other  notched  data  isusually  required  for  ASI  purposes.  When  local  plastic 
deformation  occurs,  the  local  stresses  and  strains  are  calculated  using  a  method  such  as 
Neuber'sequation.  Whentheapplied  stress  is  low,  thelocal  stresswill  si  mplybea  product  of 
the  appi  i  ed  stress  and  the  el  asti  c  stress  concentrati  on  factor.  Therefore,  theoreti  cal  I  y,  for  I  ow 
stress  and  high  cycle  fati  gue,  the  strai  n-l  ife  approach  isidentical  to  the  stress-1  ife  approach. 
FI  ence  the  strain-life  approach  may  be  considered  a  more  general  approach  for  modelling 
fatigue  damage. 

The  block  diagram  in  Figure  1  illustrates  the  process  of  the  strain-life  approach.  The  basic 
input  data  are  the  stress  concentrati  on  represent!  ng  the  geometry  of  the  speci  men,  the  cycl  i  c 
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stress-strain  relation,  the  strain-life  relation  of  the  material,  and  the  applied  load  sequence. 
The  load  sequence,  the  stress  concentration  factor,  the  stress-strain  relationship  and  the 
Neuber's  rule  are  used  to  determine  the  local  stress  and  strain  history.  Thelocal  stress  history 
is  then  processed  using  the  Rainflow  counting  method  to  form  pairs  of  peak  and  valley 
stresses  (or  stress  cycles)  that  relates  to  fatigue  damage.  The  counted  stress  cycles  are 
converted  to  strain  cycles  which  are  in  turn  converted  to  equivalent  fully-reversed  strain 
cycles,  based  on  the  selected  equivalent  strain  equation.  Finally,  the  converted  equivalent 
strai  namplitudeisusedtod  eter  mi  ne  th  e  n  u  mber  of  cy  cl  es  th  e  mater  i  alwouldhaveendured 
under  this  strain  amplitude  to  cause  failure.  The  inverse  of  this  number  gives  the  fatigue 
damage  attributable  to  this  strain  cycle.  By  summing  up  the  damages  caused  by  all  the 
Rainflow -counted  strain  cycles,  thetotal  damage  by  the  spectrumisobtained.Theinverse  of 
the  total  damage  gives  the  fati  gue  I  ife  i  n  terms  of  the  number  of  passes  of  the  spectrum.  It 
should  be  noted  that  this  approach  does  not  account  for  the  history  effect  of  the  load 
sequence.  It  uses  the  Miner  linear  summation  rule.  Thedetailsof  the  main  blocksin  Figurel 
are  discussed  in  the  foil  owing  sections. 


Figurel:  Block  diagram  illustrating  theprocess  of  strain -I  ife  approach 


From  a  conti  nuum  mechani  cs  perspective,  thestrai  n-l  ifeapproach  empi  rical  ly  relates  the  local 
elastic- plastic  behaviour  of  a  material  to  fatigue  damage.  Thus  understanding  the  elastic- 
plastic  behaviour  of  the  material  isfundamental  to  theapplicati  on  of  thestrai  n-lifeapproach. 
Given  a  cyclic  load  history,  a  methodology  isrequired  to  accurately  determinethestress and 
strain  at  the  notch  root.  FAMSFI  uses  the  Neuber's  method  which  will  bediscussed  in  the 
following  section. 

4.1  Notch  Stresses  and  Strains 

T  0  appi  y  thestrai  n-l  ife  approach,  a  means  to  accuratel  y  determi  nethestrai  n  at  the  I  ocati  on  of 
interest  is  required,  usually  the  root  of  a  notch  on  astructure.  Dueto  stress  concentrationsthat 
occu  r  at  such  I  ocati  ons,  I  ocal  stresses  may  exceed  the  yi  el  d  stress  (or  proporti  onal  I  i  mi  t)  of  the 
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material.  Under  such  conditions,  the  local  stresses  and  strains  can  no  longer  be  determined 
usi  ng  thestress  concentrati  on  factor  al  one.  A  rel  ati  onshi  p  such  as  N  euber's  equation  is  needed 
to  supplement  the  elastic  equilibrium  equations  and  the  stress-strain  relationships.  This 
simple  approach  provides  a  relatively  sound  predict!  on  of  plastic  response,  and  issuitablefor 
deal  ing  with  fatigue  loading.  Thediscussions  presented  here  are  confined  to  plane  stress  and 
uniaxial  stress  only,  which  are  the  cases  dealt  with  in  FA  MS  and  FAMSFI. 

4.1.1  N euber's  Method 

Several  methods  exist  that  may  be  used  to  estimate  the  notch  root  stress  and  strain,  but  the 
most  commonly  used  one  is  that  by  Neuber  [7].  The  strain-life  module  (FAMSFI)  in  CGAP 
uti  Uses  the  N  euber  method  for  eval  uati  ng  notch  root  stresses  and  strai  ns. 

For  elastic  deformation,  thetheoreti  cal  stress  concentrati  on  factor  can  beused  to  evaluate 

the  local  stress.  According  to  thedefiniti on  of.^:, , 


= 

^S 


(1) 


where  A5  istheremoteelasticstressrangeand  isthestressconcentration factor. Similarly 
in  the  elastic  regime  we  have 


Ae 


(2) 


where  Ae  is  the  remote  elastic  strain  range  and  is  the  strain  concentration  factor.  In  the 
el  asti  c  regi  me  both  thestress  concentrati  on  factor  and  thestrai  n  concentrati  on  factor  areequal 
to  the  theoretical  concentration  factor,  i.e. 

Neuber  noted  that  in  the  post  yield  rangethese  relations  would  no  longer  hold.  That  is,  the 
stress  concentration  factor  would  decrease  and  thestrain  concentration  factor  would 
increase  relative  to  the  elastic  stress  concentration  factor  a:,  .  N  euber  hypothesised  that  the 
elastic  stress  concentration  factor  would  beequal  to  the  geometric  mean  of  andK^.i.e., 

(3) 

Therefore  usi  ng  equation  (1)  and  (2)  equation  (3)  can  be  rewritten  to  give 


=  K^ASK^Ae  =  K^ASAe 


(4) 


As  the  remote  stress  is  assumed  to  bewithintheelasticrangewehave  A5'  =  £'Ae,  otherwise 
the  whole  cross-section  will  yield,.  Therefore 


ActAs  = 


KfAS^ 


(5) 
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where  E  is  the  modulus  of  elasticity.  This  is  Neuber's  equation  which  relates  the  remote 
el  asti  c  stress  ran  ge  to  th  e  I  ocal  el  asti  c-  p  I  asti  c  stress-strai  n  r an  ge  at  th  e  n  otch . 

Figure  2  shows  graphically  how  theNeuber  ruleisapplied  to  the  stress  strain  relationship. 
Neuber's  rulefor  a  casewhere  =l  suggests  that  the  product  of  the  applied  elasticstress 
and  strain  is  equivalent  to  the  product  of  theelastic-plasticstress  and  strain.  Geometrically, 
this  means  that  the  area  of  A  is  assumed  to  be  equal  to  the  area  of  B. 


Figure2:  A  graphical  representation  of  the  N  euber  rule  assuming  a  stress  concentration  factor  of 

4.1.2  FAMSH  Representation  of  Stress-Strain 

The  stress-strain  relationship  in  Figure  3  can  be  represented  in  various  ways.  The  Ramberg- 
Osgood  equation  is  commonly  used  to  represent  the  monotonic  stress-strain  curve,  but  can 
also  be  used  tofitthe  cyclic  stress-strain  data.  In  FAMSFI,  the  cyclic  stress-strain  curve  or 
hysteresis  stress  versus  hysteresis  strain  curve  (as  opposed  to  the  monotonic  stress-strain 
curve)  can  be  represented  in  oneof  two  ways:  either  in  tabular  or  equation  format.  I  n  tabular 
format,  the  data  may  be  cycl  i  c  stress  versus  cycl  i  c  strai  n,  or  hysteresi  s  stress  versus  hysteresi  s 
strain.  In  equation  format,  parameters  are  used  to  represent  the  cyclic  stress-strain  curve. 
Further  information  on  creating  material  definitions  can  befound  in  sect!  on  5.3.2.  More  detail 
on  the  oh  gins  of  the  FAMSFI  stress-strain  relationship  can  befound  in  the  Appendix  A. 
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Figures:  Stress-strain  reiationship  after  the  appiication  of  a  ioad  and  then  foiiowed  by  fuiiy 
reversed  ioading 


The  hysteresis  stress  versus  hysteresis  strain  curve  is  evaluated  assuming  Masing-type 
behaviour^  [8]  and  using  the  following  relations: 


Act 


(6) 

(7) 


Using theFAMSH  stress-strain  relationship  in  AppendixA  and  equations6and  7wehave 


1 


II 

2  2 

Act  f  1 
~  2E 

fAcJ 

9  ^pl 

V  ^  JJ 

>  CTpl 

(8) 

II 

Act 

'lE 

^  CTpl 

where  A^^and  As^are  the  elastic  and  plastic  components  of  the  total  strain  range  As 
respectively,  is  the  proportional  limit,  isthetotal  applied  stress  and  K'  and  «'  arethe 

cyclic  strength  coefficient  and  cyclic  strain  hardening  coefficients  respectively.  Equation  8 
describes  the  form  of  the  stabilised  hysteresis  loop  representing  the  cyclic  stress-strain 
behaviour. 

4.1.3  Fatigue  Notch  Factor 

The  previous  discussion  of  the  strain-life  method  suggests  that  the  fatigue  life  of  a  smooth 
(un-notched)  specimen  can  predict  the  fatigue  behaviour  at  the  notch  root  of  a  notched 

4  Masing-type  behaviour  assumes  that  the  hysteresis  loop  is  twice  the  cyclic  stress-strain  curve. 
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Structure.  Elastictheory  can  be  used  to  relate  the  remote  stress  to  thelocal  stress  at  the  root  of 
a  notch  using  the  stress  concentration  factor  .  In  thecase  of  fatigue  itwasoriginally  hoped 
thatthiswould  hold  trueforthefatiguestrength,  theamplitudeof  the  applied  fully-reversed 
cycl  i  c  stress  that  wi  1 1  cause  fai  I  ureafter  a  seri  es  of  repeated  cycl  es.  H  owever,  it  has  been  found 
thattheratioofthefatiguestrengthsoftheun-notched  and  notched  sped  mens  (at  the  same 
number  of  cycles)  is  not  equal  to.^,  .This  is  because  in  reality  the  notched  fatigue  strength  is 
also  affected  bythenotch  radius,  material  strength,  material  properties,  etc[9].Tocorrectthis 
behaviour  a  new  parameter,  the  fatigue  notch  factor  ,  was  introduced  to  replace.^:, . 

Thefati  gue  notch  factor  is  not  a  theoreti  cal  parameter  I  i  kethestress  concentrati  on  factor  and 
i  s  not  just  a  fu  ncti  on  of  the  geometry  and  remote  I  oad  i  ng.  Thefati  gue  notch  factor  rel  ates  the 
fatigue  strength  of  a  smooth  specimen  to  that  of  a  notched  specimen  using  the  following 
equation. 


K 


f 


Smooth  fatigue  strength 
Notched  fatigue  strength  at  equal  life 


(9) 


If  the  notched  fatigue  strength  was  just  a  function  of  the  geometry  and  the  load  then 

w  oul  d  be  equi  val  ent  to  the  el  asti  c  stress  concentrati  on  factor.  Thus  the  fati  gue  notch  factor  i  s 
essenti  al  I  y  an  empi  ri  cal  parameter  used  to  account  for  these  other  factors. 

An  empirical  relation  was  developed  by  Neuber  to  evaluate  the  fati  gue  notch  factor.  N  euber 
developed  what  is  known  as  the  N  euber  notch  factor  based  on  the  hypothesis  that  stresses 
may  be  assumed  constant  over  small  distances  [9].  This  factor  is  defined  as: 


K 


N 


=  1  + 


1  +  \la/ r 


(10) 


where  a  isthecharacteristic  length  dependant  on  thematerial  and  r  isthenotch  radius.  Some 
values  for  the  characteristic  length  a  can  be  found  in  [9].  The  Neuber  notch  factor  K,  i  sonly 
appropriate  for  unloaded  holes  or  notches  in  individual  specimens.  It  does  not  take  into 
account  any  effects  the  material  itself  might  introduce.  In  FA  M  SH ,  thefati  gue  notch  factor  is 
used  in  place  of  the  theoretical  stress  concentration  factor  to  calculate  the  local  stress  and 
strain. 


4.2  Stress  and  Strain  Related  to  Fatigue  Life 

U  nl  i  kethestress- 1  ife  approach  thestrai  n-l  ife  approach  does  not  requi  retest  data  for  a  vari  ety 
of  notch  types.  Rather,  the  strain-life  methodology  reduces  the  required  input  data  by 
providing  a  more  general  means  of  estimating  fatigue  damage  and  by  considering  prior 
knowledge  about  stress-strain  rel  ationships(i.e.  Neuber,  K,).  This  is  achieved  by  assuming 

that  the  fatigue  behaviour  at  the  notch  root  is  equivalent  to  that  of  uni-axial  test  speci  mens, 
seeFigure4.  In  other  words  thefatiguelife  of  material  at  the  root  of  a  notch  is  equivalent  to 
thefatiguelifeof  an  un-notched  specimen  of  the  same  material  subjected  to  the  same  cyclic 
straining. 
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Figure  4: 


Notch 


N  otched 
Specimen 


U  n-notched 
smooth  specimen 


Diagram  demonstrating  the  rationaie  for  using  strain-iife  data  obtained  on  smooth 
specimens  to  predict  the  fatigue  ii  fe  of  notched  specimens 


The  basic  fati  gue  I  ife  data  needed  for  the  strai  n-l  ife  approach  are  thefati gue  I  i  ves  at  vari ous 
strain  amplitudes.  These  data  are  acquired  by  conducting  tests  on  smooth  cylindrical 
specimens  subjected  to  fully-reversed  constant  amplitude  loading  under  strain  control. 
Similar  to  Basqu  in's  observation  on  stress-1  ife  data  [10],  Coffin  [11],  and  Manson  [12]  found 
thattheplastic-strain-lifedata  could  belinearised  on  a  log-log  scaleand  be  expressed  as: 

^  =  s).{2Nj.)\  (11) 


whereAf^,  is  the  plastic  strain  range,  s'j  is  the  fati  gue  ductility  coefficient,  2N j  isthenumber 

of  reversal  s  to  f  ai  I  u  re  an  d  c  i  s  th  e  d  u  cti  I  i  ty  exponent.  Combining  the  above  w  i  th  th  e  Basq  u  i  n 
equati  on  thefati  gue  I  ife  may  be  rel  ated  to  thetotal  strai  n  range  experi  enced  by  the  materi  al .  I  n 
FA  M  SH ,  thi s  rel  ati  onshi p  i s  descri  bed  i  n  the  fol  I ow i  ng  form: 


^s  Af 

—  =  — -  +  — - 
2  2  2 


2b  ^ 


Ni 


-  + 


a 


n; 


-+ 


P 


N}  N 


f 


(12) 


where  the  parameters«,  p ,  Z?'andc'are  related  to  the  fatigue  strength  and  ductility 
coeffi  ci  ents  and  exponents  vi  a: 
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E 

P^rs).  (13) 

b'^-b 

c'  --C 

For  a  particular  material  thestrain-life(£--A/ )  data  must  first  be  obtained  through  a  series  of 
tests.  These  testsinvolvetheapplicati  on  of  fully-reversed  loads  under  strain  control  to  smooth 
un-notched  specimens.  Oncethe£--A/  curveis  known  itisthen  possibleto  estimate  the  number 
of  cycles  to  failure  for  a  given  strain  history.  It  is  useful  to  recall  that  the  6--A/  diagram  deals 
with  strain  amplitude  and  as  such  the  total  strain  amplitude  is  the  sum  of  the  elastic  and 
plastic  components.  Figure  5  shows  the  total  strain  behaviour  in  terms  of  the  number  of 
reversals  to  failure  for  4340  steel  [13].  The  elastic  and  plastic  components  of  equation  12  have 
been  plotted  in  Figure  6  and  the  linear  relationship  with  fatigue  life  is  evident.  Similar 
relationships  havebeen  observed  in  many  metals. 


Figures:  Low-cyclefatiguebehaviour  of  annealed  4340  sted  showing  thetotal  strain  amplitudefor 
constant  strain  amplitude  loading[13] 
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Figures:  Low-cycle  fatigue  behaviour  of  annealed  4340  steel  showing  the  dastic  and  plastic 
components  of  strain  amplitude  for  constant  strain  amplitude  loading[13] 

4.3  FA M  SH  Representation  of  the  Strain-Life  C urve 

FAMS  and  FAMSFI  allow  a  more  flexible  representation  of  the  strain-life  curve  based  on 
equation  12.  Itisalsoworth  noting  that  FAMSFI  also  supports  tabulated  datatorepresentthe 
strain-life  curve.  While  the  codes  allow  the  user  to  implement  equation  12,  a  more  accurate 
representation  of  the  stain-life  relationship  over  the  whole  strain  range  can  be  achieved  by 
piece-wise  fitting  equation  14  to  experimental  data.  This  process  is  depicted  in  Figure  7.  The 
strain-lifedataaredivided  into  three  segments,  dominated  by  plastic  strain,  elasticstrain  or 
by  a  balanced  combination  of  elastic  and  plastic  strain.  FI  ere  the  elastic  component  of  strain 
can  be  represented  with 


As^  _  «  _  Cl  C3 

2  ~  Nf  ~  Nf  jV“  ' 

and  the  plastic  component  can  be  represented  with 

2  ~  ~  Nf  Nf  ' 


(14) 


(15) 


In  Figure7ELET and  PLET arestrain  rangeschosentodefinethetransitions between  regions 
1,  2  and  3.  Depending  on  the  values  of  E  LET  and  PLET,  different  combinations  of  equations 
areused  to  represent  theenti  re  strain-life  curve,  asdetailed  in  Table!  The  constants  used  in 
Table  1  and  Figure  7  can  be  found  by  curve-fitting  equation  12  to  experimental  data.  See 
A ppendix  B  for  an  exampi e  of  determi  ni  ng  the  i  nput  parameters  for  the  FA  M  SFI  strai  n-l  ife 
equation. 
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Reversals  to  Failure  IN^ 


Figure  7;  Strain-lifecurverepresentation  in  FAM  S  and  FA  M  SH 


T abiel:  Parameters  used  in  FA  M  SH  to  defined  each  region  of  thestrain-iife  curve  in  Figure  7 


Total  strain 

Region  1 

Region  2 

Region  3 

PLET=ELET=0 

f(Ci,Bi,C2,Gi) 

N/A 

N/A 

ELET>PLET=0 

^(Ci,Bi,C2,Gi) 

^(C3,B2,C2,Gi) 

N/A 

PLET>ELET=0 

^(Ci,Bi,C2,Gi) 

^(Ci,Bi,C4,G2) 

N/A 

ELET>PLET>0 

^'(Ci,Bi,C2,Gi) 

^(C3,B2,C2,Gi) 

^(C3,B2,C4,G2) 

PLET>ELET>0 

^(Ci,Bi,C2,Gi) 

^(Ci,Bi,C4,G2) 

^(C3,B2,C4,G2) 

4.4  Rainf low  Counting 

Rainflow  counting  is  an  internal  function  that  is  conducted  by  the  FAMSH  code.  It  is 
performed  on  the  local  stress  history  derived  from  the  applied  load  spectrum  using  the 
N  euber  rul  e.  The  I  oad  history  is  i  nput  i  nto  FA  M  SFI  as  a  seri  es  of  turni  ng  poi  nts.  Thestrai  n-l  ife 
algorithm  used  in  FAMSFI  requires  that  this  spectrum  be  Rainflow  counted.  Rainflow 
counting  hasaphysical  basisand  thatiswhen  metals  are  subjected  to  repeated  loading  the 
stress-strain  response  forms  hysteresis  I  oops  I  ike  those  in  Figures.  Iftheloading  remainsin 
the  elastic  regime  then  the  hysteresis  I  oops  remain  closed  such  as  that  produced  by  thecycle 
G  to  FI .  If  the  loads  applied  are  large  enough  to  cause  plastic  yielding  then  open  hysteresis 
loops  are  formed.  It  is  useful  to  note  that  if  the  loading  does  not  cause  plasticity,  then  any 
hysteresis  I  oops  formed  degeneratetoastraightline.  Any  actual  loop,  with  non-zero  enclosed 
area,  i  nd  i  cates  pi  asti  city.  Thearea  contai  ned  by  the  hysteresi  s  I  oops  represents  the  energy  I  ost 
during  that  cycle.  In  the  case  of  a  variable  amplitude  I  oad  sequence,  hysteresis  loops  can  be 
d raw n  such  that  onef i  nds  smal  1 1  oops  contai  ned  w i  thi  n  I  arger  I  oops.  The  Rai  nfl  ow  counti  ng 
scheme  developed  by  M  atsuishi  and  Endo  [14]  may  be  used  to  extract  the  peak  and  valley  of 
all  the  closed  hysteresi  si  oops  from  a  complex  spectrum.  It  is  based  on  the  observed  memory 
effect  in  many  metals  where  a  smaller  load  cycle  appears  to  contribute  an  interruption  to  a 
larger  load  cycle.  In  Figure  8  the  load  cycle  B-C  serves  as  a  smal  I  interruption  to  the  larger 
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load  cycleA-D.  At  point  C  as  the  load  i  ncreases  back  to  B  the  stress-strain  curve  resumes  its 
ori  gi  nal  path  from  A  to  D  i  n  effect  rememberi  ng  the  previ  ous  I  oad  appi  i  cati  on  from  A  to  B. 

TheRainflow  counting  method  can  be  envisaged  by  rotatingthestrain-timehistorysothatthe 
verti  cal  axi  s  i  s  the  ti  me  and  the  hori  zontal  axi  s  the  strai  n  ampi  i  tu d e.  The  strai  n  h i  story  can  be 
vi  ewed  as  a  seri  es  of  si  opi  ng  roofs  that  overl  ap  one  another  at  d  i  fferent  poi  nts.  The  Rai  nfl  ow 
cycles  are  then  defined  by  the  manner  in  which  rain  is  allowed  to  drip  down  the  roofs.  A 
number  of  rules  are  used  to  define  the  behaviour  of  therein. 

The  process  of  Rai  nfl  ow  counting  isasfollows[14] : 

•  Allow  therein  to  start  from  the  largest  peak  (or  valley)  and  run  down  the  roof  till  it 
falls  off  completely. 

•  Start  the  next  d  ropi  et  of  water  from  the  val  I  ey  (or  peak)  at  the  poi  nt  w  here  the  d  ropi  et 
fell  from  the  roof  and  allow  it  to  run  in  the  opposite  direct! on  till  it  falls  off  the  roof. 
This  constitutes  one  Rai  nfl  ow  counted  cycle. 

•  Next  begin  again  from  the  first  load  reversal  that  has  not  been  fully  touched  by  a 
droplet.  Allow  the  droplet  to  run  till  it  falls  off  and  terminate  this  path  when  either  it 
falls  off  completely  or  hits  a  previous  path  created  by  a  previous  droplet. 

•  F rom  thi  s  reversal ,  al  I  ow  a  d ropi  et  to  f I  ow  i  n  the  opposi  te  d  i  recti  on  and  term!  nate  i  n  a 
similar  manner.  This  now  becomes  the  next  Rainflow  counted  cycle.  Continue  the 
process  until  all  the  cycles  within  theload  spectra  havebeen  counted. 

In  Figure  8  there  are  4  full  Rainflow  cyclesA-D,  B-C,  E-F  and  G-FI  which  correspond  to  the 
4 closed  stress-strain  hysteresis  loops.  For  further  information  refer  to  Fuchs,  Fatemi  and 
Stephens  [15]  or  Ellyin  [16]. 
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Figured:  Pictorial  Representation  oftheRainflow  Counting  M  ethod 

4.5  Mean  Stress  Effects  and  Equivalent  Strain  Equations 

For  fully  reversed  strain  cycles  (i.e.  R  =  mini  mum  stress/  maximum  stress  =-l)  the  strain 
range  can  bedirectly  used  to  obtain  the  corresponding  fatigue  life  from  the  6--A/  curve.  Fora 
vari  abl  e  I  oad  hi  story,  the  strai  n  cycl  es  are  usual  I  y  not  fu  1 1  y  reversed ,  hence  i  t  i  s  i  mportant  to 
account  for  the  effect  of  (non-zero)  mean  stresses  on  fati  gue  I  ife.  A  n  equati  on  that  converts  a 
general  strain  cycle  where  -I  i  nto  an  equ  i  valent  cycl  e^^^  where  i?  =  -l  isknownas 

an  equivalent  strain  equation,  and  several  models  have  been  developed  over  the  years. 
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I  n  general  Dow  I  i  ng  [  17, 18]  provi  des  some  gui  dance  to  the  use  of  equ  i  val  ent  strai  n  equati  ons. 

H  e  conci  udes  the  fol  I  owi  ng: 

•  The  M  orrow  equation  is  reasonably  accurate  in  most  cases  when  the  true  fracture 
strength  is  known. 

•  The  Smith  Watson  and  Topper  equation  providesgood  resultsin  most  cases  and  for 
aluminium  alloys  produces  somewhat  more  accurate  resultswhen  compared  to  the 
Morrow  equation. 

•  In  cases  where  there  is  enough  data  to  determine  the  exponent  min  the  Walker 
equation,  the  Walker  equation  gives  superior  results. 

•  For  high  strength  aluminium  alloys  the  exponent  min  the  Walker  equation  is 
approximately  0.5  reducing  the  Walker  equation  back  to  the  Smith  Watson  Topper 
equati  on.  For  I  ower  strength  al  umi  ni  um  al  I  oys  the  exponent  is  hi  gher. 

FAMSFI  implements  the  fol  I  owing  eight  equivalent  strain  equations: 

i)  Modified  Morrow  equation  [19] 


s 


eq 


I _  ^ mean 


(Tp 


(16) 


where  is  the  mean  stress  and  isthetruefracturestrengthstheengineeringstressat 
the  beginning  of  fracture  during  a  tensile  test. 

ii)  Loopin  equation  [20] 

Both  the  oh  gins  of  the  Loopin  and  modified  loopin  equati  ons  are  described  in  [20]. 


^eq  ^a 


^eq 


(  \P 

2  _l_  ^  mean 


V 


'a  J 
\ 


2 _  ^mean 


C^mean  ^  0 


C^mean  <  0 


a  J 


(17) 


where  cj^  istheamplitudeof  thelocal  stress  and  p  isa  material  parameter. 


5  True  fracture  strength  is  defined  as  the  load  at  failure  divided  by  the  actual  cross-sectional  area  at 
failure 
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iii)  Modified  Loopin  equation  [20] 


^eq 


^  a 

2  _l_  mean 


^eq 


V  J 

f  \-aip 

j _  ^mean 
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O-mean  ^  0 


*^mean  ^  ® 


where  «  and  p  are  material  parameters, 

iv)  Walker  equation  [21] 


^eq 


\  1— A« 


(18) 


(19) 


where  is  the  peak  stress  and  m  is  a  material  parameter. 

v)  Smith  Watson  Topper  (SWT)  equation  [22] 

P  _  P  ‘^max 

^eq 

\  E 

The  Smith  Watson  and  Topper  mean  stress  equation  is  a  general  formulation  that  has  been 
show  n  to  correl  ate  w el  I  w  i  th  the  fati  gu e  exper i  ments  of  al  u  mi  n i  u  m  al  I  oys  [ 23] . 

vi)  Modified  F-18 equation  [22,  24] 

The  report  by  Ghi  del  I  a  et.  al .  [24]  provi  des  a  good  background  to  thedevel  opment  of  the  F-18 
equivalent  strain  equation.  An  important  point  to  note  is  that  the  equation  was  developed 
using  data  that  did  not  include  negative  mean  stresses. 


O-max  ^  0 


C^max  <  0 


(20) 
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^eq^ 
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frr  ^ 

+ 

^mean 

V 

1  rE  j 

\-R 


O-mean  ^  0 


<^mean  <  0 


where  f  is  a  material  parameter. 


(21) 


4.6  Fatigue  Damage  Accumulation 

The  fatigue  damage  caused  by  one  cycle  is  defined  by  the  reciprocal  of  the  number  of 
repetitions  of  that  cycle  needed  to  cause  failure.  The  loss  of  energy  per  hysteresis  loop  is 
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ad  d  i  ti  ve  w  h  i  ch  I  eads  to  the  assu  mpti  on  that  fati  gue  d  amage  i  s  al  so  ad  d  i  ti  ve  (i  .e.  the  Pal  mgren- 
Minerrule[25]).Thislineard  amage  accu  mu  I  ati  on  concept  w  as  f  i  rst  proposed  by  Pal  mgren  i  n 
1923  [26]  and  Miner  provided  the  mathematics  to  describeltin  1945  [25]. 

Consequently,  the  fatigue  damage  produced  byn  cycles  at  one  strain  level  Is  defined  as 

£»  =  —  (22) 

N 


where 

n  Isthenumber  of  applied  cyclesatstrain  rangelevel  As^-, 

N  Isthetotal  numberof  cycles  to  cause  failure  at  the  strain  rangelevel  read  off 

the  stral  n-l  Ife  curve  for  the  given  stress  range  I  evel . 

If  the  spectrum  consists  of  repeated  cycles  at  one  constant  strain  amplitude,  then  failure  Is 
assumed  to  occur  when  n  =  N  and  thus  the  damage  D  =  l 

For  a  varl  abl  e  ampi  I  tude  I  oad  history  contal  nl  ng  more  than  onestral  n  range,  fatl  gue  damage 
can  be  calculated  asthesumofthedamageattrlbutabletoeach  Individual  strain  rangelevel 
within  thespectrum.  That  Is 


sa=i: 


N, 


Failure  Is  assumed  to  occur  when 


z 


D  >  1 


(23) 


(24) 


Th  I  s  1 1  near  su  mmatl  on  of  damage  caused  by  d  I  fferent  stral  n  ranges  I  s  know  n  as  the  Pal  mgren- 
Mlner  rule[25].To  apply  this  summation  rule,  thespectrum  must  first  be  analysed  to  Identify 
the  closed  hysteresi  s  I  oops  using  a  scheme  such  as  Rainflow  counting.  Each  hysteresis  loop  Is 
converted  to  an  equivalent  strain  cycle  with  zero  mean  stress  using  an  appropriate  equivalent 
strain  equation.  Following  this  step,  thePalmgren-Mlner  rule  equation  23  can  be  applied  to 
obtain  the  total  damage  resulting  from  the  spectrum.  Once  the  damage  calculated  using 
equation  23  exceeds  1,  failure  of  the  structure  Is  assumed  to  occur.  In  the  FAMSFI 
I  mpl  ementatl  on,  the  I  nverseof  thetotal  damage  I  n  equatl  on  23 1  s  used  to  cal  cul  ate  the  number 
of  times  (passes)  the  spectrum  Is  applied,  such  that  thetotal  damage  Is  1.0.  This  number  Is 
then  used  to  determine  the  total  fatigue  life  In  flight  hours  using  the  flight  hour's  one 
spectrum  represents. 

Presently  FAMSFI  only  ImplementsthePalmgren-Mlner  ruletocalculatefatiguedamagebut 
It  must  be  remembered  that  there  areother  techniques  of  fatigue  damage  calculation.  Since 
1945  a  pi  ethora  of  methods  have  been  proposed  to  eval  uatefatl  guedamage,  mal  niy  to  account 
for  the  observed  effect  of  the  load  sequence  on  thetotal  I  Ife  of  structures  [27]. 
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5.  Using  FAMSH  in  CGAP 

5.1  Starting  FAMSH 

WhenCGAPislau  nched  from  the  M  i  crosoft  W  i  nd  ows  Start  men  u ,  the  user  i  s  presented  w  i  th 
the  screen  as  shown  in  Figure  9.  By  defauitthefirstmoduledisplayed  is  the  CGAP  Solver  for 
crack  growth  analysis. 


Figure  9:  D  efault  Start  up  mode  for  CGAP 


To  use  the  strain-life  model  the  user  must  first  switch  to  theFAMSFI  module  in  the  CGAP 
GUI.  This  can  be  achieved  in  two  ways: 

Either  by  pressing  the  "Configure"  button  located  inthetoptool  bar  or  by  selecting  "tools  ^ 
configuration..."  on  the  menu  bar.  The  following  "Settings"  dialog  box  should  appear 
Figure  10. 
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Figure  10:  TheCGAP  sditings  dialog  box 

In  this  dialog  box  click  the  drop  down  box  titled  "Solver  Module  Name"  and  select 
"FAMSH".Then  click  "OK". 

FAMSH  should  automatically  linktothedatabasefileon  installation,  however  the  database 
can  be  selected  or  reselected  at  a  later  stage  by  using  the  "Open  DatabaseFile"  option  under 
the  "File"  menu. 

TheFAMSH  moduleshould  now  beloaded  and  you  are  ready  to  begin.  The  moduleshould 
I  ook  I  i  ke  Fi  gure  11.  The  words  "Runni  ng  FA  M  SH  "  shoul  d  appear  i  n  the  bottom  I  eft  panel  to 
indi cate  that  the  FA  MSFI  module  has  loaded  properly. 


Figurell:  FAMSH  module  loaded  in  CGAP 
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5.2  Data  Inputs  and  Outputs 

TheFilenamainputfilecan  be  created  manually  and  used  as  an  input,  butCGAP  will  create 
this  fi  I e  for  the  user  w hen  the  user  i  nput  the!  r  data  i  nto  the  GU I . 

The  input  files  used  to  specify  a  problem  prior  to  analysis  are  as  follows: 

•  Filename.!  nput 

•  Filenamaspectra 

•  Ftn07  (Required  only  for  DSA^ formats) 

•  Ftnl5  (Optionally,  the  single  value  contained  in  this  file  may  be  entered  in 
Filenamainput,  on  the  last  line  with  the  keyword  "  Maximum  Delta  Stress:” 
followed  by  a  value  such  as  60000) 

The  output  files  that  are  created  by  theCGAP  FAMSFI  moduleare: 

•  Filenamainput 

•  Filenamapasses 

•  Filenamatruncspc 

•  Filenamadsout 

•  Filenamadmo 

•  Flmatl.dat 

•  Flmat2.dat 

See  sect!  on  5.8  for  more  detail. 

There  are  two  ways  to  enter  input  data: 

•  Open  an  existing  input  file,  ag.,  from  previous  analysis; 

•  Enter  data  directly  through  CGAP  GUI. 

In  either  case  the  main  input  file  has  the  extension  ".input"  and  it  contains  all  the  control 
information  for  the  anal  ysis.Thespectrumfilehastheextensi  on  ".spectra"  and  itcontainsthe 
load  spectrum,  defined  by  turning  points.  Thedamagesourcecodefileisentitled  "ftn07"  and 
it  contains  the  codes  for  different  missionscomprisingthespectrum,  see  reference  [2].  These 
codes  are  called  "Damage  Source  Assignments"  or  DSAs^.  Optionally,  a  file  "ftnl5" 
contai  ni  ng  a  si  ngl  e  number  represent!  ng  the  maxi  mum  stress  range  i  n  the  spectrum  can  be 
used.  The  ".input"  file  can  be  created  with  CGAP  by  entering  data  into  CGAP  directly  and 
pressi  ng  the "  B u i  i  d  "  button .  The spectru m  f i  I  e  and  the  d amage  sou rce  cod e f i  I  e  are  normal  I  y 
generated  by  upstream  software  but  can  also  be  created  using  a  text  editor.  The  material 
specifications  are  contained  in  either  the  "fimatl.dat"  (tabulated  data)  or  the  "fimat2.dat" 
(equati  on  speci  fi  cati  on)  f  i  I  es.  These  fi  I  es  are  generated  automati  call  y  w  hen  the  user  sel  ects  a 


6  Damage  Source  Assignments  were  developed  for  the  P-3C  SLAP.  See  reference  [2]  for  more 
information. 

1  The  Damage  Source  Assignment  was  a  modification  made  by  Lockheed  Martin  to  examine  the 
proportions  of  damage  associated  with  different  in  flight  manoeuvres. 
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material  from  the  material  database  in  the  CG  A  P  GUl.Thesefilescan  be  used  as  inputs  using 
theoriginal  FAMSH  executable. Thesinglenumber  contained  in  "ftnl5"  isused  to  generate 
thedamagedistribution  table,  and  does  not  affect  thecal  culated  fatigue  I  ife.Thestress  levels 
are  categorised  i  nto  20  bi  ns  usi  ng  this  val  ue.  By  default  CG  A  P  chooses  60000  if  no  number  has 
been  specified. 

5.2.1  Units 

To  maintain  consistency  and  backward  compatibility  with  older  versions  of  the  FA  MS  and 
FAMSFI  codes  it  is  recommended  that  the  foil  owing  units  be  used  when  entering  data  into 
CGAP. 

Table2:  Recommended  input  units  for  CGAP 


Input  parameter 

Units 

Young's  Modulus 

ksi 

Fracture  Strength 

ksi 

Proportional  Limit 

ksi 

Ultimate  Strength 

ksi 

PTQI 

Vksi 

A 1 1  stresses  i  n  the  stress-strai  nor  stress-1  i  fe  tabi  e 

ksi 

Residual  Stress  (Resd) 

ksi 

MAXRNG 

psi 

Spectrum  file 

psi 

5.2.2  Importing  an  Analysis  File 

Importing  data  can  be  accomplished  using  the  "Open"  button  on  the  top  tool  bar  or  by 
clicking  "File^Open". 

The  i  nput  fi  I  e  i  s  a  text  docu  ment  w  i  th  the  extensi  on "  .input" .  TabI  e  3  gi  ves  an  exampi  e  of  the 
input  file  and  an  associated  list  of  the  source  code  variables. 
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Tables:  An  eompleof  the  input  data  file  use  in  CGAP 


Line 

No. 

Example.input 

FAMSH 

Variables 

Line  description 

1 

generic  test  seq 

JTITLE 

Input  title 

2 

1 

IPTl 

M  aterlal  data  output  information 

3 

12 

IPT2 

Analysis  Output  1  nformati  on 

4 

1 

IPT3 

Loads  Data  Format  (ASCI  1  or  Binary) 

5 

0 

IPT4 

Damage  Source  1  ndicator 

6 

0 

IPT5 

N  umber  of  Passes 

7 

1 

KN 

Notch  Factor 

8 

D7075-T651 

MATERIAL 

Material  Identification 

9 

1 

MATE 

Material  Data  Type 

10 

1 

UNITCV 

Unit  Conversion 

11 

0 

RS 

Residual  Stress 

12 

2 

NMLT 

N  u  mber  of  Stress  M  u  1  ti  p  1  i  ers 

13 

2.5  5 

SREF 

Stress  Multipliers 

14 

le-007 

DCUT ( 9 ) 

Damage  Truncation  Level 

15 

2 

Equivalent  strain  equation 

16 

Flight  Hours :  1 

Spectrum  file  total  flight  hours* 

17 

Maximum  Delta 
Stress : 60000 

Originally  in  aseparatefile"ftnl5.dat"=** 

*  Here,  Flight  Hours :  and  Maximum  Delta  Stress :  are  keywords  and  must  be  entered  as  shown. 

**  Lines l-14mustbeentered  inthesequenceshown.  Lines  15-17 may  beentered  inanyorder,  orany  linemay 


be  omitted.  When  omitted,  the  equivalent  strain  equation  defaults  to  2,  theflight  hours  default  to  land  the 
maximum  delta  stress  defaults  to  60000. 

InputTitJe 

On  the fi  rst  I i ne  any  i nformati on  can  be  i nput  up  to  a  I ength  of  80  characters. 

Material  Data  Output  I  nformati  on 

This  line  contains  an  integer  number  that  controls  the  type  of  material  data  written  to  the 
output  data  fi  I  e.  The  fol  I  ow i  ng  opti  ons  are  aval  I  abl  e: 

0  No  materi al  output  -  N o  stress strai n  data  is  output 

1  Cycl  i  c  Strai  n  Stress  -  A  strai  n  vs  stress  tabi  e  i  s  pri  nted  to  the  output  fi  I  e 

2  H  ysteresi  s  Stress  Strain  -  A  stress*strai  n  vs  stress  tabI  e  i  s  pri  nted  to  the 

output  file 


A  nalysi  s  0  utput  I  nformati  on 

Another  integer  number  is  used  to  control  the  amount  of  information  that  is  output  to  the 
output  data  fi  I  e.  The  fol  I  ow  i  ng  opti  ons  are  aval  I  abl  e: 


-11 

or 

-1 

Basic 

10 

or 

0 

Detailed 

11 

or 

1 

More  Detailed 

12 

or 

2 

Most  Detailed 

Ifthesecond  number  above  is  used  then  no  Tpassesfile  is  produced.  For  more  information 
see  section  5.5.3. 
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Loads  Data  Format  (A  SC  1 1  or  Binary) 

This  line  contains  an  integer  number  that  tells  the  program  whether  the  spectrum  file  is  in 
Binary  or  ASCI  I  format 
0  Bi  nary  format 

1  ASCII  format 

D  amage  Source  I  ndicator 

This  is  an  integer  number  that  tells  the  program  whether  the  Damage  Source  Assignment 
(DSA)  code  is  used.  The  DSA  is  a  nine  digit  number  that  is  used  in  the  spectrum  file  and  is 
used  inthesubsequent  anal  ysisand  output[2].  For  further  information  and  an  exampleofthe 
application  of  damage  source  assignments  refer  to  [28]. 

0  DSA  is  used 

1  DSA  is  not  used 

N  umber  of  Passes 

An  integer  number  is  used  to  switch  between  one  pass  and  two  passes.  A  residual  stress  is 
calculated  after  one  pass  and  is  used  in  the  calculation  of  fatigue  life  during  the  second  pass. 
0  ne  pass  w  i  1 1  i  gnore  any  resi  d  ual  stress  cal  cu  I  ated  after  the  fi  rst  pass. 

0  Two  passes  required 

1  One  pass  required 

N  otch  Factor 

The  notch  factor  istheNeuber  notch  factor  Kj^  and  not  the  stress  concentration  factor  a:,  .  It 
should  be  noted  that  scaling  factors  are  applied  to  this  number  using  a  different  input 
parameter  to  be  discussed  later. 

Material  Identification 

This  variable  is  a  character  string  of  length  25.  In  earlier  versions  the  character  string  was 
limited  to  9  characters.  This  variable  specifies  the  material  to  be  used  and  reflects  the  exact 
name  used  i  n  the  material  database.  If  the  namedoes  not  exist  i  n  the  database  a  warni  ng  will 
be  displayed  by  theCGAP  GUI. 

Material  Data  Type 

MATF  isan  integer  that  indicates  whether  material  data  is  specified  in  table  format  or  as  an 
equati  on .  I  f "  1"  i  s  sel  ected  then  the  program  expects  the  mater i  al  data  to  be  i  n  a  tabi  e  I  ook-u  p 
format  in  afilewith  the  name  "flmatl.dat".  If  "2"  isselected  then  the  program  expects  the 
material  data  to  be  in  the  form  of  an  equation  specified  in  the  "flmat2.dat"  file.Thiswill  be 
automatically  selected  by  CGAP  using  the  information  in  thedatabase. 

1  Table  look-up 

2  Equation 

Unit  Conversion 

Thislinecontainsareal  numberwhich  isusedtoconverttheunitsofthematerial  datafile.To 
convert  from  KSI  to  M  Pa  use  the  conversi  on  factor  of  6.895.  It  does  not  affect  the  I  oads  i  n  the 
spectrum  file. 


23 


DSrO-TR-2392 


Residual  Stress 

Thisline  contains  the  magnitudeof  the  residual  stress  needed  in  the  analysis.  If  no  residual 
stress  is  to  be  used,  then  enter  0.0.  The  units  of  the  resi  dual  stress  must  be  consistent  with  the 
units  of  the  proportional  limit,  that  isksi. 

N  umber  of  M  ultipliers 

This  integer  parameter  is  used  to  specify  the  number  of  Stress  multipliers  to  use.  The 
multi  pi  i  ers  are  specifi  ed  on  the  next  I  i  ne  of  the  i  nput  data  fi  I  e. 

Stress  Multipliers 

An  array  is  used  tostorea  maximum  of  10  real  numbers.  Each  number  in  the  array  specifies  a 
mu  I  ti  pi  i  cati  on  factor  to  use  i  n  conj  uncti  on  w  i  th  the  notch  factor.  The  notch  factor  i  s  mu  I  ti  pi  i  ed 
by  one  of  the  stress  mu  I  ti  pi  i  ers  w  h  i  ch  i  n  tu  rn  i  s  used  to  seal  e  the  appi  i  ed  spectru  m.  A  nal  ysi  s 
occurs  for  each  of  the  stress  multipliers.  If  the  notch  factor  is  given  a  value  of  1  then  stress 
multipliers  can  be  used  to  specify  a  series  of  different  notch  stress  concentration  factors. 

Truncation  Level 

Thislinecontainsareal  number  that  is  used  to  denotetheuser  specified  truncation  level.  This 
truncation  level  truncates  the  spectrum  based  on  the  associated  damage  calculated  for  each 
cycle.  Thenumber  used  hereisadamagelevel  below  which  cycleswith  damageless  than  this 
value  are  ignored. 

Equivalent  Strain  Equation 

An  integer  number  is  used  to  select  the  equivalent  strain  equation  to  use  in  the  analysis. 
Several  equivalent  strain  equations  are  availablefor  the  user  and  are  as  follows.  Section  4.5 
describes  these  equations  in  greater  detail. 

2  M  odi  fi  ed  M  arrow  equ  ati  on 

The  next  options  required  extra  information  to  be  placed  in  the  ".input"  file.  The  extra 
information  isentered  on  the  nextiinefollowingtheequivalent  strain  equation  option  inthe 
".input"  file. 

3  Loopin 

The  beta  parameter  in  the  Loopin  equation  is  required  on  the  next  line 
underneath  3 

4  M  odi fied Loopin 

Two  parameters  used  in  the  Modified  Loopin  equation,  beta  and  alpha,  are 
required  on  the  next  line  underneath  4 

5  1/1/  a/ker 

A  parameter  "m'  used  in  the  Walker  equation  is  required  on  the  next  line 
underneath  5 

6  5  nnith-Watson -Topper 

A  parameter "  m'  used  i  n  the  Smith- Watson-T  opper  equati  on  is  requi  red  on  the 
next  I  i  ne  underneath  6 
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7  ModifiedF-18 

Two  parameters  used  in  the  Modified  F-18  equation,  beta  and  gamma,  are 
required  on  the  next  line  underneath  7 

Spectrum  FileTotal  Flight  Hours 

Thi  s  n  u  mber  sped  f  i  es  the  n  u  mber  of  f  I  i  ght  hou  rs  the  spectru  m  represents.  The  text  "Flight 
Hours;"  mu  St  be  entered  exactly  as  shown  followed  by  the  number  of  flight  hours. 

M  aximum  D elta  Stress 

The  I  argest  change  I  n  stress  I  n  the  spectru  m  i  s  used  I  n  the  cal  cu  I  ati  on  of  damage  d  ensi  ti  es.  The 
text  "Maximum  Delta  Stress : "  must  be  entered  exactly  as  shown  fol  I  owed  by  a  number. 
Thevalue  here  is  not  critical  tothefinal  life  estimates  produced.  Thefuncti  on  of  this  variable 
influencesthe  bin  sizes  used  in  determining  the  percentage  of  damage  caused  by  different 
loads.  This  information  is  recorded  in  thefilenamadsoutouput. 

There  are  three  tabs  that  are  used  to  enter  data.  They  are  the  "Material",  "Load"  and  "Case 
C  ontrol "  tabs. 

5.3  Material  tab 

The  materi  al  tab  (see  Fi  gure  12)  al  I  ows  the  user  to  create  a  new  materi  al ,  or  I  oad  or  update  an 
existing  material  in  the  database.  The  material  properties,  including  Young's  modulus,  the 
proportional  limit,  the  fracture  strength,  the  ultimate  strength,  the  stress-strain  relationship 
and  thestrai  n-l  ife  rel  ati  onshi  p  are  aval  I  abl  eto  al  I  the  cases  usi  ng  thesame  database.  This  tab 
also  all  ows  the  user  to  select  the  strain-1  ife  equation,  equivalent  strain  equation  and  to  enter 
related  parameters.  These  data  are  not  stored  inthedatabase;they  aresaved  in  the  input  file, 
and  hence  only  availablefor  the  current  case  being  analysed. 
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FigureU:  Thematerialstab  used  in  theFAM  5H  CGAP  module 


5.3.1  Loading  an  Existing  Material 

Materials  from  the  database  can  beselected  using  the  Material  Databasedropdown  menu,  see 
Figure  13. 


Material  |  Load  |  Case  Contro 


Material  Database 


-d 

2014-T6 

2024T351 

2024-T851 

7050-7651 -Eqn 

7050-T7351 

ksi 

7050-T74 

|7075-T651 

ksi 

Figure  13:  The  material  drop  down  menu 


5.3.2  Creating  a  New  Material  Database  Entry 


To  create  a  new  material,  enter  a  new  name  or  modify  an  existing  material  ID  in  the 
Material  Database  dropdown  box.Thiswill  enabletwodropdown  menus:  "Stress- Strain 
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Data"  and  "Data  Type".  The  user  may  now  select  either  the  tabular  input  or  the  equation 
input  option  from  the  "Stress-Strain  Data"  dropdown  menu,  see  Figure  14. 


Stress-Strain  Data 

T  abular 

Equation  | 

T  abular  ■ 

Figure  14:  Stress-Strain  Data  dropdown  menu 

The  "Data  Type"  dropdown  menu  (see  Figure  15)  can  be  used  to  specify  the  form  of  the 
stress-strain  relationship.  This  option  only  affects  data  entered  in  tabular  form;  it  does  not 
affect  the  equation  format.  The  equation  parameters  always  define  the  cyclic  stress-strain 
rel  ati  onsh  i  p .  Th  ree "  D  ata  Ty pe"  opti  ons  are  avai  I  abl  e 

•  Flys  sts  V  hys  stn  x  hys  sts  (Flysteresis  stress  versus  hysteresis  strain  x  hysteresis 
stress) 

•  Cyc  stress  v  eye  strai  n  (Cycl i  c  stress  versus  cycl  i c  strai  n) 

•  FI  ys  stress  v  hys  strai  n  (hysteresi s  stress  versus  hysteresi s  strai  n) 


Figure  15:  D  ata  Type  dropdown  menu 

Data  can  be  entered  directly  into  the  grids  provided,  as  shown  in  Figurel6.  The  number  of 
data  isspecified  for  both  the  stress-strain  curves  and  the  strain-life  curves. 


Material  Properties 

Number  of  Data  Stress-Strain  Data 

I  I  □ 


Cyclic  Strain 

Cyclic  Stres!  ^ 

1 

0.00587 

60.5 

2 

0.00612 

63 

3 

0.00636 

65.5 

4 

0.0066 

68 

5 

0.00685 

70.5 

<  > 

Data  Type  [ 


N  umber  of  D  ata  S  train-Life  D  ata 

P  I  □ 


-  I  EqStnEq|  ^ 


Figure  16:  M  ateriai  tabuiar  input  in  the  FA  M  SH  CGAP  moduie 
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Data  entered  i  n  the  tabi  e  can  be  del  eted  by  ri  ght  cl  i  cki  ng  the  specific  data  you  w  ish  to  del  ete 
and  selecting  "Delete",  see  Figure  17. 


Insert  Above 

Insert  Below 

tyclic  Stres!  ^ 

120 

Delete 

Iso 

— 1 - 

34 

35 

3030 

340 

3580 

350 

36 

4200 

360 

rt 

1<  > 

Data  Type  [" 


Figure  17:  D  deling  data  in  themateriai  data  tabies 


M  aterial  data  can  also  beinput  in  theform  of  equations.  Selecting  "Equation"  from  the  Stress- 
Strai  n  Data  d  ropdow  n  menu  w  i  1 1  al  I  ow  the  user  to  enter  the  rel  evant  constants  that  specify  the 
material  equations  used  in  FAMSFI .  An  example  of  the  material  property  specification  is 
shown  in  Figure  18. 


Material  Properties 


Stress-Strain  Data 


1 - □ 

C1 

0.01039 

C2  1 

0.17734 

C3 

0.02204 

C3  [ 

5316 

PTQI  fT 

16519e-0  1/ksi 

Q  1 

3.5 

Data  Type  [ 


Strain-Life  Data 


1 - 

B1 

0.07363 

82  r 

0.1645 

G1 

0.57856 

G2  r 

2.1493 

ELET 

0.005747 

PLET  P 

0.0103 

I  EqStnEq]  ^ 


Figure  18:  M  ateriai  equation  input  in  theFAMSH  CGAP  moduie 


The  stress-strai n  equation  used  in  FAMSFI  (equation  8)  is  presented  in  section  4.1.2.  In 
FAMSFI  the  input  parameters  used  arePTQI  and  Q  and  rel  atetothecyclic  strength  coefficient 
K'  and  cyclic  strain  hardening  coefficients  n'  via  the  foil  owing  equations: 

PTQI^^  (25) 

A 

e=^  (26) 

n 
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Refer  to  section  4.3  for  a  detailed  description  of  the  constants  used  in  FAMSH  to  define  the 
strain-life  equation.  Appendix  B  explains  how  to  determine  these  parameters  from 
experimental  data. 

Oncethedata  has  been  entered,  the  new  material  specification  must  be  saved  tothedatabase. 
To  do  this  press  the  save  button  (see  Figure  19)  next  to  the  Material  N  ame  i  nput  box. 


Figure  19:  C lick  the  save  button  to  save  the  material  to  the  database 

To  remove  the  material  from  thedatabase  use  thebutton  marked  with  acrossnexttothesave 
button  in  Figure  19. 

5.3.3  Eq  Stn  Eq  (Equivalent  Strain  Equation) 

The  equ  i  val  ent  strai  n  equ  ati  on  i  s  al  so  sel  ected  u  nd  er  the  materi  al  s  tab.  The  avai  I  abl  e  opti  ons 
are  detailed  in  section  4.5. 

5.3.4  Strain-Life  Equation 

This  menu  is  for  information  only  and  is  not  selectable  by  the  user.  Its  function  is  to  indicate 
which  strain-life  equation  is  being  used  in  conjunction  with  the  equivalent  strain  equation. 
When  this  box  either  contains  the  word  "equation"  or  remains  empty  then  the  default  stain 
life  equation,  equation  12  is  used.  If  it  shows  "Tabular",  then  the  tabular  strain-life  data  is 
used. 

5.3.5  UnitCv  (Unit Conversion) 

The  vari able  UnitCv  is  used  to  convert  the  units  of  the  material  properties.  For  example,  to 
convert  from  the  Imperial  unitksi  to  the  SI  unitM  Pa,  a  value  of  6.895  should  beused.  ltdoes 
not  affect  the  loads  in  the  spectrum  file. 

5.4  Load  tab 

The  I  oad  tab  (see  Fi  gure  20)  al  I  ows  the  user  to  sel  ect  the  spectrum  fi  I  e  and  to  enter  a  number 
of  load-related  parameters. 
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Figure20:  The  load  tab  used  in  the  FA  M  5H  CGAP  module 

5.4.1  Spectrum  File... 

Thespectrum  file  is  loaded  using  the  "Spectrum  File... "  button.  By  defauitthespectrumfile 
has  the  extension  ".spectra". 

5.4.2  Spectrum 

The  user  must  i  nd  i  cate  if  the  spectrum  fi  I  e  contai  ns  DSA  i  nformati  on;  w  hether  i  t  i  s  i  n  bi  nary 
or  text  format.  Thi  s  i s  done  by  sel  ecti  ng  the  appropri  ate  opti  on  u nder  the  I  oad  tab  usi  ng  the 
Spectrum  dropdown  menu.  There  are  four  options  aval  I  able  and  shown  in  Figure  21. 

Spectrum  [Binary  Damage  Source  Assignme  -^1  I” 


Binary  Damage  Source  Assignment 


Binary  No  Damage  Source  Assignmer 
ASCII  Damage  Source  Assignment 
ASCII  No  Damage  Source  Assignmen 

Figure21:  Spectrum  dropdown  menu 
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A  fifth  option  is  available  "DBI/  SST"  formats.  This  can  be  selected  by  choosing  "ASCII 
D  amageSourceA  ssignment"  and  checking  the  "D  BI/SST  Format"  box  next  to  the  "Spectrum" 
dropdown  menu.  For  further  information  on  DBI/  SST  format  refer  to  [6]. 

5.4.3  MAXRNG  (Maximum  Delta  Stress) 

The  number  entered  here  is  the  maximum  applied  stress  range,  which  is  used  to  calculate 
damage  distribution.  The  maximum  stress  range  can  be  calculated  by  subtracting  the 
mi  n  i  mu  m  spectru  m  stress  from  the  maxi  mu  m  spectru  m  stress,  but  cu  rrenti  y  a  val  u  e  of  60000 
isused  by  default.  Itshould  denoted  that  thisvaluedoes  not  affect  thecal  culatedfatiguelife. 
Thi  s  val  ue  w  i  1 1  be  cal  cu  I  ated  automati  cal  I  y  i  n  the  next  versi  on  of  FA  M  SFI . 

5.4.4  Scale 

This  was  originally  referred  to  as /Cn  in  the  FAMScode,  but  has  been  renamed  to  Scale  in 
CGAP  to  better  reflect  its  functionality.  This  value  is  essentially  a  scaling  factor  to  the 
spectrum.  The  actual  applied  stress  is  the  product  of  this  value,  the  level  ofloading  and  Kn. 

5.4.5  Resd  (Residual  Stress) 

Residual  stress  can  be  introduced  using  the  "Resd"  input  under  the  load  tab.  The  residual 
stress  should  use  the  same  units  as  the  proportional  limit. 

5.4.6  Dmg  (Damage  cut  off) 

This  is  the  damage  cut  off  threshold  DCUT.  This  specifies  the  truncation  level.  That  is  it 
truncates  loads  in  the  spectrum  with  damage  less  than  DCUT. 

5.4.7  Flours  per  sequence 

This  value  is  used  to  convert  the  life  in  cycles  to  flight  hours.  This  number  appears  in  the 
".input"  file  after  the  key  phrase  "Flight  Hours:" 


^TheDBI/  SSI  format  wasdevdoped  to  supporttheP-SCsavicelifeassessment  program  and  was  used 
to  examine  the  damage  resulting  from  different  flight  conditions. 
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5.5  Case  Control  tab 

The  Case  Control  tab  (see  Fi  gu re  22)  al  I  ows  the  user  to  sel  ect  appropri  ate  output  opti  ons. 


Figure22:  The  Case  Control  tab  used  in  the  FA  M  SH  CGAP  module 


5.5.1  Problenn  Title 

Under  this  tab  the  user  can  add  a  title  to  the  output  usi  ng  the  "Problem  Title"  input  box. 

5.5.2  N  0  sref_vs_passes/ sref_vs_passes 

The  user  can  also  chooseto  create  a  ".passes"  file  by  toggling  between  "Nosref_vs_  passes" 
and  "sref_vs_  passes".  See  section  5.8for  more  information. 

5.5.3  Analysis  Output  Options 

There  arefour  options  under  the  "Analysis  Output  Options"  dropdown  menu.  The  user  has 
the  opti  on  to  choose  from "  Basi  c" , "  D  etai  I  ed  " , "  M  ore  D  etai  I  ed  "  and "  M  ost  D  etai  I  ed  " .  These 
opti  ons  control  the  amou  nt  of  i  nformati  on  w  ri  tten  to  the "  .output"  fi  I  e  d  u  ri  ng  anal  ysi  s.  Basi  c 
level  information  reports  the  total  damage  calculated  and  the  total  life  as  well  as  other  low 
level  model  information.  Detailed  information  additionally  reports  on  the  damage  each  block 
in  the  spectrum  produces.  More  Detailed  output  also  includes  the  Rainflow  cycles  in  the 
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spectru  m  and  the  notch  stress  and  strai  n  cal  cu  I  ated  for  each  of  these  cycl  es.  The  most  detai  I  ed 
option  also  gives  information  about  the  notch  stress  and  strain,  the  alternating  strain,  mean 
stress,  equivalent  stress,  failure  cycles  and  the  damage  for  each  cycle. 

5.5.4  Material  Data  Output 

Under  the  "Material  Data  Output"  dropdown  menu  the  user  can  choose  what  material  data 
to  output  to  the  output  fi  I  e.  The  user  has  the  choi  ce  of "  N  o  materi  al  output" ,  "C  yd  i  c  Stress- 
Strain"  and  "HysteresisStress-Strain".The"CydicStress-Strain"  option  outputs  a  table  of 
datadescri  bing  the  cyclic  stress-strain  curve.  The  "Hysteresis  Stress-Strain"  option  outputs  a 
tabi  e  of  d  ata  th  at  d  escr  i  bes  th  e  hysteresi  s  str  ess-strai  n  cu  rve. 

5.5.5  Passes 

The"Passes"  dropdown  menu  provides  the  user  the  option  to  use  one  or  two  passesinthe 
analysis  process. 

5.6  Spectrum  file  formats 

There  are essenti ally  fivespectrumfileformatsavailableand  they  are: 

•  Bl ock  format  with  Desi gnated 

•  Damage  Source  Assignment(DSA) 

•  Bl  ock  format  w  ithout  Damage  Source  Assi  gnment  (DSA ) 

•  Cycleformat  with  Damage  Source  Assignment  (DSA) 

•  Cycleformat  without  Damage  Source  Assignment  (DSA) 

•  DBI/ SST  input  format 

The  format  of  the  spectru  m  f  i  I  es  i  s  i  mportant  i  f  the  anal  ysi  s  i  s  to  r  u  n  correcti  y .  The  f  ol  I  ow  i  ng 
sections  describe  in  detail  theformatto  use  for  each  of  the  spectrum  file  formats. 
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5.6.1  Block  format  with  D  amage  Sou  rce  Assign  merit  (DSA) 

The  layout  and  format  of  the  spectrum  file  that  is  required  when  using  the  Damage  Source 
Assignment  and  Block  type  inputs  is  presented  inTable4. 


T able  4:  T he  format  of  a  FA  M  SFI  input  file  using  D  SA  and  Block  formats  (fixed  format) 


ExampleDSABIock.input 


2 

1 

BLOCK 

The 
number 
of  flights 

T 

Arbitrary 
number 
that  is 
not  used 

T 

Spectrum 

format 

1  38 

FLIGHT  NUMBER  &  MISSION  TYPE 

Y 

Flight 

number 

Y 

Mission 

Type 

General  descrip^n  not  used  in 
the  analysis 

2 

1.000  TURNING  POINT  PAIRS  &  FLIGHT  DURATION 

^ 

Y 

Number 
of  lines  in 

Y 

Flight 

Duration 

General  description  not  used  in 
the  analysis 

the  flight 


38001022 


-13408 


"Y" 


DSA  Maximum  Minimum  Number 
TP  stress  TP  stress  of  TPPs 


38003019 


6891 


1087 


Each  number  should  be 
^  contained  within  the 
allocated  spaces 
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5.6.2  Block  format  without  Damage  Source  Assignment  (DSA) 

Thelayout  and  format  of  thespectrumfilethatisrequired  when  using  Block  type  inputs  and 
no  DSA  is  presented  inTableS.  Note  that  the  arbitrary  number  seen  in  the  DSA  block  format 
in  the  first  line  should  be  removed.  Also  note  that  the  maximum  and  minimum  stresses 
should  be  expressed  using  integer  numbers  and  not  real  numbers. 


T able  5:  T he  format  of  a  FA  M  SH  input  file  using  the  Block  format  with  no  D  5A 


Example  no  D  SA  Block-input 


1 


The 
number 
of  flights 


BLOCK 

_ _ ^ 

V 

Spectrum 

format 


Example  Flight  Spectra 

Y 

TitI e  of  the  spectrum 


2  1.000  TURNING  POINT  PAIRS  &  FLIGHT  DURATION 

Number  Flight  General  description  not  used  in 

ofTPPs  Duration  the  analysis 

in  the 


flight 

1  -13408  4 

'  Y  T 

Maximum  Minimum  Number 
TP  stress  TP  stress  ofTPPs 


6891  1087  2 

T  Y  y 

10  spaces  10  spaces  10  spaces 
Integer  Integer  Integer 


Each  number  should  be 
^  contained  within  the  10 
spaces 
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5.6.3  Cycle  format  with  DamageSourceAssignment  (DSA) 

The  layout  and  format  of  the  spectrum  file  that  is  required  when  using  the  Damage  Source 
Assignment  and  Cycletype  inputs  is  presented  inTableG. 


T able  6:  T he  format  of  a  FA  M  5H  input  file  using  D  5A  and  Cycle  formats 


ExampleDSACycle.input 


1 


The 
number 
of  flights 

1 


Flight 

number 


Arbitrary 
number 
that  is 
not  used 
38 


- Y - 

Mission 

Type 


CYCLE 

T 

Spectrum 

format 


FLIGHT  NUMBER  &  MISSION  TYPE 


General  description  not  used  in 
the  analysis 


114 


Number 
ofTPPs 
in  the 


1.000  TURNING  POINT  PAIRS  &  FLIGHT  DURATION 


Y 

Flight 

Duration 


General  description  not  used  in 
the  analysis 


flight 

38001022  1.00  38001022  -13408.00 


Y 

Y 

Y 

Y 

DSA 

Maximum 

DSA 

Minimum 

TP  stress 

TP  stress 

38002023 

J 

-1773.00 

V  ) 

38002023 

V  ) 

-9068.00 

V  i 

\  ^  \ 

L  Y  1 

12  spaces 

10  spaces 

10  spaces 

10  spaces 

Integer 

Real 

Integer 

Real 

The  user  has  the  opti  on 
^  to  enter  another  TPP  on 
the  same  line 


I  Each  number  should  be 
^  contained  within  the 
allocated  spaces 


One  Space  One  Space  One  Space 
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5.6.4  Cycle  format  without  Damage  Source  Assign  merit  (DSA) 

The  layout  and  format  of  thespectrumfilethat  is  required  when  using  Cycletype  inputs  and 
no  DSA  is  presented  inTable?.  Note  that  the  arbitrary  number  seen  in  the  DSA  block  format 
in  the  first  line  should  be  removed. 


T  able  7;  T  he  format  of  a  FA  M  SH  input  file  using  the  Cycle  format  with  no  DSA 


ExamplenoD  SACycle.input 


1 


The 
number 
of  flights 


CYCLE 

_ _ _ _ ^ 

V 

Spectrum 

format 


Example  Flight  Spectra 

Y 

Title  of  the  spectrum 


Y 

Number 
ofTPPs 
in  the 
flight 


114 


1.000  TURNING  POINT  PAIRS  &  FLIGHT  DURATION 


~Y 

Flight 

Duration 


General  description  not  used  in 
the  analysis 


1.00  -13408.00 


Maximum 
TP  stress 


Minimum 
TP  stress 


-1773.00 

Y 

Maximum 
TP  stress 


-9068.00 

j 

Y 

Minimum 
TP  stress 


The  user  has  the  opti  on 
^  toenter  upto4TPPson 
the  same  line 


6891.00  1087.00 


10  spaces  10  spaces 
Real  Real 


15104.00 


Y 

10  spaces 
Real 


5740.00 

Y  ' 

10  spaces 
Real 


> 


J 


Each  number  should  be 
contained  within  the  10 
spaces 
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5.6.5  DBI/ SST  in  put  format 

The  layout  and  format  of  the  spectrum  filethat  is  required  when  using  the  DBI/  SST  type 
inputs  is  presented  in  Tables. 


T able  8:  T he  format  of  a  FA  M  5H  input  file  using  D  BI/55T  formats 


ExampleDBICycle.input 


1 


The 
number 
of  flights 

1 


Flight 

number 


Arbitrary 
number 
that  is 
not  used 
38 


- Y - 

Mission 

Type 


CYCLE 

T 

Spectrum 

format 


FLIGHT  NUMBER  &  MISSION  TYPE 


General  description  not  used  in 
the  analysis 


114  1.000  TURNING  POINT  PAIRS  &  FLIGHT  DURATION 

- '' - . - 

Number  Flight  General  description  not  used  in 

ofTPPs  Duration  the  analysis 

in  the 
flight 


0005001091  869.62  0005001091 


DSA  Maximum  DSA 
TP  stress 


295.63  0005004072  3884.06  0005004072  -1068.51 


Minimum  dSA  Maximum  DSA  Minimum 
TP  stress  TP  stress  TP  stress 


0005004072  3884.06  0005004072  -1068.51  0005004072  4819.06  0005004072  -1690.94 


Y  i 

Y  ^ 1 

j 

i  Y  i 

j 

i  Y  i 

k  Y  , 

^  j 

k  Y 

10  spaces 
Integer 

10  spaces 
Real 

10  spaces 
Integer 

10  spaces 
Real 

10  spaces 
Integer 

10  spaces 
Real 

10  spaces 
Integer 

10  spaces 
Real 

One  Space  One  SpaceOne  Space  One  Space  One  Space  One  Space  One  Space 
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5.7  Analysis  Using  FAMSH 

Performi  ng  the  anal  ysi  s  usi  ng  CGA  P  i  s  rather  si  mpl  e.  0  nee  the  rel  evant  data  has  been  input 
intoCGAP,  all  the  user  needs  to  do  isto  press  the  "Build"  button  to  build  the  case  and  then 
press  "Run"  to  executethe  analysis,  see  Figure  23. 


Rie  Edit  View  Commands  Tools  Help 

i  -  H  X 

New  : .  Save  Remove 

ib  Mp 

Import  Export  Save  Out 

Build  Run  Stop 

\ 

0 

Configure 

Figure23:  CGAP  menu  bar 

A  cti  ve  anal  ysi  s  i  s  i  nd  i  cated  by  the  progress  bar  at  the  bottom  of  the  w  i  ndow ,  see  Fi  gu  re  24. 


Progress  llllllllllllllllllli 

Figure24:  Progress  bar 

When  the  analysis  is  complete,  a  message  will  be  displayed  in  the  messages  pane  at  the 
bottom  of  thewindow. 

5.8  FAMSH  output 

Several  outputs  are  created  byFAMSFI.A  bri  ef  description  of  each  of  the  output  filesisgiven 
i  n  TabI  e  9.  ExampI  es  of  the  output  are  provi  ded  i  n  A  ppend i  x  C . 


T able  9:  FA  M  5H  output  files 


Output  File 

Description 

Example 

filenamedmo 

The  .dmo  file  contains  the  damage  matrix.  The  damage  matrix 
provides  the  user  with  information  relating  to  the  distribution  of 
damage  based  on  the  mission  types.  This  information  is  only 
available  when  the  DSA  format  is  used. 

Table  23 

filenamedsout 

This  file  records  the  damage  density  table  calculated  using  the 
maxi  mu  m  change  i  n  stress  defi  ned  i  n  the  ftnl5  f  i  1  e  or  i  nput  f  i  1  e.  The 
.dsoutfileisequivalenttotheftnlBfileproduced  by  FAMS. 

Table  24 

filenamepasses 

The  first  line  of  this  file  contains  the  path  of  the  input  file.  The 
fol  1  ow  i  ng  1  i  nes  contai  n  the  Kn  values  specified  intheinputfileusing 
the  source  variable  SREF  along  with  the  evaluated  fatigue  life  in 
fl  i  ght  hou  rs  to  fai  1  u  re. 

Table  25 

filenamatruncspc 

Thisfilecontainsthedamage based  truncated  spectrum. Thisfileis 
equivalent  to  theftnll  file  produced  by  FAMS. 

Table  26 

filenameoutput 

This  file  contains  all  the  output  produced  by  FAMSFI.  The 
i  nformati  on  i  n  the  output  f  i  1  e  i  s  control  1  ed  by  the  case  control  tab  i  n 
CGAP. 

Table  27 
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6.  Material  Database 


CGAP  incorporates  a  new  material  database  that  links  seamlessly  with  the  CGAP 
environment  for  crack  growth  and  fatigue  analysis.  The  database  allows  the  user  to  select, 
retrieve  and  store  material  properties  to  perform  analyses.  The  use  of  an  integrated  material 
database,  in  comparison  to  individual  material  datafiles,  provides  greater  consistency.  This 
thereby  minimises  errors  associated  with  incorrect  material  input  files  and  improves 
traceability  of  the  analysis  results.  The  database  also  provides  support  to  quality  control  of 
outputs  of  the  analyses  for  certification  and  validation  purposes. 

6.1  Materials  Data  Structure 

A  diagram  of  the  CGAP  materials  data  structure  is  provided  in  Figure  25.  The  fields  in  the 
four  data  tabi  es  pertai  ni  ng  to  crack  i  niti  ati  on  properti es  are  descri  bed  i  n  TabI e  10. 


Figure  25:  CGAP  M  aterial  D  ata  Structure  Diagram 
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TableName 

TabI  e  D  escri  pd  on 

CGAP  Database 

Field 

FAMSH  V»i^le 

Description 

Material_CI 

Contains  unique 
materiai  identifier  as 
weil  as  common 
materiai  properties. 
Links  with  aii  the 
other  crack  initiation 
parameter  tabies. 

id 

MTITLE(1:9) 

M  ateri  ai  speci  f i  cati  on 

name 

MTITLE(10:45) 

Materiai  description 

MatFormat 

MATF 

The  materiai  format 

M  atF  ormat  =  1 T  abu  i  ar 

MatFormat  =2  Equation 

E 

E 

Young's  Modulus 

Sproportionai 

PLMT 

proportional  limit 

Sultimate 

FULT 

ultimate  strength 

NPTH 

Sfracture 

fracture  strength 

N_pts_CycHys 

SF 

the  number  of  cyclic/  hysteresis  data  pairs 

N_pts_EN 

NPTS 

the  number  of  materiai  strain-life  data  pairs 

CycHysType 

IM 

a  f  1  ag  that  speci  f i  es  the  type  of  cycl  i  cj  hysteresi  s  d  ata 

IM  =1  -  cyclicstrain  versusstressdata 

1 M  =  2  -  hysteresi  s  strai  n  versu  s  stress  d  ata 

1 M  =  3  -  hysteresi  s  stress-strai  n  versu  s  stress  d  ata 

factor 

FCTR 

the  factor  required  to  convert  data  from  ksi  into  psi 

Material_CycHys_CI 

Contains  crack 
initiation  data  for 
cyclic  /  hysteresis 
stress-strain  curves 
represented  by 
discrete  data  points 
(from  ffmatl.dat file). 
Links  with 

Material_CI  and 

M  ateriai_StrainLife_CI 
tabies. 

id 

MTITLE(1:9) 

M  ateri  al  speci  f  i  cati  on 

CycHysData 

HYSl(NPTH) 

contains  either  cyclic  strain,  hysteresis  strain  or  hysteresis  stress- 
strain  data. 

StressData 

HYS2(NPTH) 

contai  ns  stress  data. 
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TableName 

TabI  e  D  escri  pti  on 

CGAP  Database 

Field 

FAMSH  V»i^le 

Description 

M  ater  i  al  _Strai  n  L  i  f e_C  1 

Contains  crack 
initiation  data  for 
material  strain-life 
curves  represented  by 
discrete  data  points 
(from  ffmatl.dat file). 
Links  with  the 
Material_CycHys_CI 
tableand  Material_CI 
tables. 

id 

MTITLE(1:9) 

M  ateri  al  sped  f  i  cati  on 

StrainData 

SNLl(NPTS) 

contains  material  strain  amplitude  data. 

LifeData 

SN  L2(N  PTS) 

contains  material  lifedata  in  number  of  cycles. 

M  ateri  al_Equati  on_C  1 

Contains  crack 
initiation  data  for  the 
stress-strain 
relationship  and  the 
strain-life  curves 
represented  by 
equations  (from 
ffmatzcfetfile).  Links 
with  the  M ateri al_CI 
table. 

id 

MTITLE(1:9) 

M  ateri  al  sped  f  i  cati  on 

PTQLQ 

PTQLQ 

Parameters  used  to  describe  the  cyclic  stress-strain  curve(see 
section  4.1.2  for  more  information) 

Cl,  C2,  C3,  C4 

C1,C2,C3,C4 

Parameters  related  to  the  fatigue  strength  and  ductility 
coefficients  and  exponents.  See  equation  12  and  equation  13.  Also 
seeAppendix  B  for  a  suggested  method  to  evaluate  these 
coefficients. 

Bl,  B2,  Gl,  G2 

Bl,  B2,  Gl,  G2 

Exponents  related  to  the  fatigue  strength  and  ductility  exponents 
via  equation  13.  Also  seeAppendix  B  for  a  suggested  method  to 
eval  uate  these  exponents. 

PLET,  ELET 

PLET,  ELET 

are  strains  chosen  to  define  the  transitions  between  region  1,  2 
and  3 
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6.2  FAMSH  interaction  with  the  Material  Database 

The  material  properties  stored  in  the  database  are  used  by  FAMSH  indirectly  through  one  of 
two  material  files,  flmatl.dat  and  flmat2.dat,  automatically  generated  by  the  CGAP  GUI.  A 
sample  of  each  of  these  files  is  provided  in  Appendix  D.  Both  files  contain  basic  material 
property  data  such  astheYoung's  M  odulus,  yield  and  ultimate  strengths.  They  also  contain 
material  cyclic/ hysteresis  curves  and  strain-life  curves.  They  differ  in  that  flmatl.dat 
represents  the  material  curves  via  discrete  data  points,  while  flmat2.dat  uses  up  to  three 
equations  to  represent  the  cyclic/  hysteresis  and  strain-life  curves.  A  "flag"  variable(MATF) 
within  the  FAMSH  input  file  tells  the  program  which  material  data  file  to  utilise.  The 
parameters  and  format  of  flmatl.dat  and  flmat2.dat  are  presented  in  Table  11  and  Table  13 
respectively.  The  definitions  of  the  variables  contained  in  flmatl.dat  and  flmat2.dat  are 
presented  i  n  TabI  e  12  and  T abl  e  14,  respecti  vel  y . 

Table  11:  Format  oftheflmatl.dat  file 


Line 

File  Parameters 

1 

P  MTITLE 

2 

E  PLMT  FULT  SF 

3 

NPTH  FCTR  IM 

4 

HYSl(l),  HYS2(1)HYS1(2),  HYS2(2). . .  HYSl(NPTH),  HYS2(NPTH) 

5 

NPTS 

6 

SNL(l),  SNL2(1)  SNL1(2),  SNL2(2) . . .  SN  L1(N  PTS),  SN  L2(N  PTS) 

Tablel2:  Flleparameter  descrlptlons(flmatl.dat) 


File  Parameter 

Description  of  the  file  parameter 

P 

a  single  character  used  by  FAMSH  todistinguishbetweenacomment(represented  bya";") 
or  a  material  header  line  (represented  by  a  ">") 

MTITLE 

material  title 

E 

Young's  Modulus 

PLMT 

proportional  limit 

FULT 

ultimate  strength 

SF 

fracture  strength 

NPTH 

the  number  of  cyclic/  hysteresis  data  pairs 

FCTR 

the  factor  required  to  convert  data  from  ksi  into  psi 

IM 

a  f  1  ag  that  speci  f i  es  the  type  of  cycl  i  cj  hysteresi  s  d  ata 

IM  =1- cyclicstrain  versusstressdata 

1 M  =2  -  hysteresis  strain  versus  stress  data 

1 M  =  3  -  hysteresi  s  stress-strai  n  versu  s  stress  d  ata 

HYSl(NPTH) 

array  with  either  cyclic  strain,  hysteresis  strain  or  hysteresis  stress-strain  data.  Forms  data 
pairs  with  HYS2(NPTH). 

HYS2(NPTH) 

array  with  stress  data.  Forms  data  pairs  with  HYS1(N  PTH). 

NPTS 

the  number  of  material  strain-life  data  pairs 

SNLl(NPTS) 

array  with  material  strain  amplitude  data.  Forms  data  pairs  with  SN  L2(N  PTS). 

SNL2(NPTS) 

array  with  material  life  data  in  number  of  cycles.  Forms  data  pairs  with  SN  L1(N  PTS). 
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T able  13:  Format  of  the  flmat2.dat  file 


Line 

File  Parameters 

1 

P 

MTITLE 

2 

E 

PLMT 

FULT 

SF 

PTQI 

Q 

3 

Cl 

C2 

C3 

C4 

4 

B1 

B2 

G1 

G2 

PLET 

ELET 

T able  14:  File  parameter  descriptions(flmat2.dat) 


File  Parameter 

Description  of  the  file  parameter 

P 

a  singlecharacter  used  by  FAMSFI  to  distinguish  between  acomment  (represented  by  a";") 
or  a  material  header  line  (represented  by  a  ">") 

MTITLE 

material  title 

E 

Young's  Modulus 

PLMT 

proportional  limit 

FULT 

ultimate  strength 

SF 

fracture  strength 

PTQLQ 

Parameters  used  to  describe  the  cyclic  stress-strain  curve  (see  section  4.1.2  for  more 
information) 

Cl,  C2,  C3,  C4 

Parameters  related  to  the  fatigue  strength  and  ductility  coefficients  and  exponents.  See 
equation  12  and  equation  13 

Bl,  B2,G1,  G2 

Exponents  related  to  the  fatigue  strength  and  ductility  exponents  via  equation  13 

PLET,  ELET 

are  strains  chosen  to  define  the  transitions  between  region  1,  2  and  3 

N  ote  that  FA  M  SH  has  I  i  mi  ts  on  the  si  ze  of  the  arrays  associ  ated  w  i  th  the  fimatl .  dat  f  i  i  e.  A  rrays 
FIYSI  and  HYS2  can  contain  no  more  than  75  eiements  whiie  arrays  SNLl  and  SNL2  can 
contain  no  more  than  65  eiements.  CGAP  hasa  iimit  of  50  eiements. 


7.  Example 

This  section  uses  an  exampie  to  iiiustrate  the  use  of  the  FAMSFI  code.  It  is  helpful  to 
remember  that  the  FAMSFI  module  operates  in  a  very  similar  manner  to  the  FAMS  code. 
Thus,  wecan  expect  similar  outputs  to  beproduced.  Foradetailed  description  with  regard  to 
thestrain-lifetheory  on  which  FAMSFI  isbased,  refer  to  section  4.  The  example  presented  in 
this  section  will  becovered  in  four  sections  which  will  explaintheProblem,  Inputs,  Execution 
and  Results. 

7.1  Problem 

Consider  a  problem  in  which  the  same  spectrum  is  applied  at  two  locations  on  a  structure 
manufactured  from  7050-T7351  aluminium.  The  locations  on  the  structure  will  be  denoted 
I  ocati  on  1  and  I  ocati  on  2.  Both  I  ocati  on  1  and  I  ocati  on  2  havegeometri  c  features  that  result  i  n  a 
stress  concentration.  Assumethatweknow  the  Neuber  stress  concentration  factors(.^^ )  in 
these  locations  are  2.5  and  5  respectively. 

For  the  purposes  of  this  exampi  e  consi  der  a  spectrum  composed  of  thefol  I  owi  ng  sub-bl  ocks 
of  constant  amplitude  loading,  shown  in  Table  15. 
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T able  15:  Example  spectrum  loading 


Constant 
amplitude  sub¬ 
block  number 

DSA 

R  ratio 

^max(PSi) 

Number  of  Cycles 

1 

22 

-1 

12000 

5 

2 

19 

0.5 

12000 

20 

3 

5 

0.8 

10000 

100 

4 

19 

0.5 

12000 

20 

The  aim  of  the  analysis  is  to  determine  the  estimated  time  to  failure  at  location  1  and 
location  2. 

7.2  Inputs 

Wewill  call  thisexample"Gen_seq"  and  name  all  our  input  files  accordingly.  To  perform  the 
analysis  we  need  to  creates  input  files.  Thesefiles  are  the  ".input",  ".spectra"  and  "ftn07" 
files. 

Thefirst  step  is  to  create  the  spectrum  file.  Forthisexamplethespectrumfilewill  bein  ASCII 
format,  in  blocks  and  using  DamageSource  Assignments,  which  can  be  created  using  any  text 
editor.  Section  5.6.1  presents  the  description  of  this  format  for  the  spectrum  file. 

T abl  e  16  shows  the  compi  ete spectrum  fi  I  e.  Li  ne  1  contai  ns  three  numbers.  Thefi  rst  number, 
"  1",  i  ndi  cates  onefi  i  ght  has  been  defi  ned  i  n  thespectrum  fi  I  e.  The  second  number  is  a  dummy 
variablewhich  isnotused  in  the  analysis.  Enter  "1"  as  the  second  number.  Finally  thekey 
word  "BLOCK"  should  beentered  to  indicate  thatthespectrumfileisdefi  ned  in  block  format. 
Line2contai  nstwo  numbers,  thefi  i  ght  number  and  the  mi ssi  on  type.  Thefi  i  ght  number  is  not 
used  in  the  analysis  but  is  used  to  indicate  which  flight  the  analysis  refers  to.  Enter  "101"  for 
thefi  ight  number.  The  mission  type  identifies  themi  ssi  on  usi  ng  a  number.  In  this  examplewe 
are  identifying  the  mission  as  mission  27  so  enter  "27"  as  the  second  number  on  line  2.  Any 
textthatisadded  following  thesenumbers  will  not  be  used  byFAMSFI  in  the  analysis. 

L  i  ne  3  al  so  contai  ns  tw  o  nu  mbers;  the  fi  rst  nu  mber  i  s  the  total  nu  mber  of  I  i  nes  bel  ow  L  i  ne  3 
that  defi  nethespectrum,  and  thesecond  number  defi  nes  thefi  i  ght  durati  on.  Si  nee  four  I  i  nes 
are  needed  in  thisexampleto  describe  the  flight,  weenter  "4"  for  thefirst  number.  For  this 
problem  let  usfurther  assume  that  thespectrum  represents  lOOfI  ight  hours  so  enter  "100"  for 
the  second  number.  Although  100  is  entered  here  it  is  not  used  in  the  analysis.  The  actual 
flight  hours  represented  by  thespectrum  is  input  separately  in  the  .input file. 

The  next  four  I  i  nes  contai  n  the  i  nformati  on  shown  i  n  TabI  e  15.  For  thi  s  probi  em  name  thefi  I  e 

"Gen_seq. spectra"  and  save  it. 
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T able  16:  T he  spectrum  file 


G  en spec.spectra 

1 

1 

BLOCK 

101 

27 

FLIGHT 

NUMBER  &  MISSION  TYPE 

4 

100.0 

CYCLE 

COUNT  &  FLIGHT  DURATION 

27001022 

12000 

-12000.0 

5 

27002019 

12000 

6000.0 

20 

27003005 

10000 

8000.0 

1000 

27004019 

12000 

6000.0 

20. 

As  this  problem  uses  DSA's  it  requires  an  extra  file  to  describe  the  DSA  codes.  This  file  is 
called  "ftn07"  and  is  also  a  text  file.  For  this  problem,  create  a  file  with  the  following 
information  in  Table  17  and  save  it  with  the  name  "ftn07". 


Table  17:  Theftn07file 


Ftn07 

'MISSIONS'  27 

' DAMAGE  SOURCES ' 

3 

005  'R=0.8 

019  'R=0.5 

022  'R=-l 

To  create  the  ".input"  file  you  can  either  create  it  using  a  text  editor  or  using  CGAP.To  use 
CGAP  open  CGAP  and  enter  the  FA  MSFI  module  as  described  in  section  5.1.  By  defauitthe 
material  stab  should  bedisplayed.  Under  the  materialsdatabasedropdown  box  select  7050- 
T7351  as  shown  in  Figure  26. 


Material 


Material  Database 


Load 


Case  Centre 


7050-T7351 


3  hIkJ 


201 4TG  .V 

2014-T6a 

'^2014-T6b 

HRJPREC 

2014TGC 

2014-TGk 

_ 

2024-T351 

ksi 

2024-T851 

ksi 

7050-7G51-Eqn 

I7050-T7351  ■ 

ksi 

I7050-T74  1 

Figure 26:  M  aterials  tab  M  aterials  D  atabase  drop  down  menu 


The  materi  al  s  tab  w  i  1 1 1  ook  I  i  ke  Fi  gu  re  27  w  i  th  the  tabu  I  ar  materi  al  d  escri  pti  on  acti  ve.  Both  the 
Cyclic  Stress- Strain  Curveand  Strain-LifeCurvefigures should  now  depopulated. 
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Mftieiial  I  Load  |  Case  Control  |  Output 


Material  Name  |ALALLOY(NORTHR)PBEC 
Static  Properties 

Young's  Modius  |1COOO  ksi 

Fractue  Strength  90  Lsi 

Proportional  UmK  40  ksi 

Ulimate  Strength  ksi 


Ur^CorMt  [t" 


Matenal  Properties 

Number  of  Data  Stiess-Sliain  Data  Number  of  Data  StraivUfe  Data 

I -  r — 3  I  r— — 3 


C3 


Cyck  Stran 

Cydk  Stres:  ^ 

Life 

IsbainArrp 

1 

0 

0 

1 

10 

0.085 

2 

0.64 

80 

2 

20 

O.OSl 

3 

1,09 

100 

3 

30 

0.039 

4 

1.43 

no 

4 

40 

0.0325 

5 

1.82 

118 

5 

60 

0.025 

< 

> 

6 

86 

0.02 

Data  Type  |  -I  Eq  Stn  Eq  |  Modeled  Motrow 


Cyclic  Stress-Strain  Curve 


0  3.68  7.36  11  04  14  72  18.4 

2xStress  x  2xStrain 


Graph  | 


Stiai>Life  Equation  I  3 


IE-3  .  1  p  1 

1E1  ie2  163  1E4  165  1E6  167 


Life  in  Cycles 


Figure 27:  M  aterials  tab  with  thetabuiar  materiai  description  active 


Now  click  on  the"Eq  Stn  Eq"  dropdown  menu.  Select  "Modified  Morrow"  for  equivalent 
strain  equation.  Enter  lin  the"UnitCv"  box. 

Now  select  the  load  tab.  Enter  the  following  information  (refer  to  section  5.4  for  more 
information  on  what  each  input  is). 

•  In  theSpectrum  dropdown  box  select  "ASCII  Damage  Source  Assignment". 

•  Under  Hours  per  Sequence  enter  100. 

•  Under  "MAXRNG"  enter  30000. 

•  Under  "Scale"  enter  1 

•  Under  "Resd"  enter  0 

•  Under  "D mg"  enter  le-7 

•  I  nthesectioncontainingthetablewiththetitle"Str«s  Levels",  therearetwo  columns 
titled  "No."  and  "Kn".  Inthefirstrow  under  "Kn"  enter  the  of  location  1  which  is 
2.5.  In  the  next  row  enter  the  of  location  2  which  is  5.  You  will  note  that  the 
nu  mber  next  to  the  box  ti  tl  ed "  Stress  Level  s"  has  i  ncremented  to  2  i  nd  i  cati  ng  there  are 
two  cases  to  analyse. 

•  Whenselectingthespectrumfiletheusercanchoosetoselectthenameand  path  using 
the  "Spectrum  file"  button,  see  Figure  28.  Ifthespectrumfilehasthesamenameand 
exists  in  the  same  directory  as  the  input  file  then  a  path  and  name  do  not  need  to  be 
specified.  In  this  case  just  delete  the  default  name  "cstamp"  and  leave  it  blank; 
otherwise,  click  on  "Spectrum  File"  tobrowsethefilesystemand  select  the  spectrum 
file. 
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Figure28:  Spectrum  file  input 

Now  select  the  Case  Control  tab  and  enter  the  foil  owing. 

•  Select  "serf_vs_ passes"  and  under  the  "N  umber  of  Passes  to  Analyse"  dropdown 
menu  select  "First  Pass". 

•  Under  the  "Material  Data  Output"  dropdown  menu  select  "Cyclic  Stress-Strain". 

Thisoption  will  output  cyclic  stress-strain  data  to  the  output  file. 

•  Leave  all  other  options  as  default. 

After  you  haveentered  all  the  above  data  the  Case  Control  tab  should  look  like  Figure  29. 


Material 


Load 


Case  Control 


Output 


Case  Database 


U  a  XI 


Problem  Title 

I  Crack  Initiation  Life  Analysis 

Input  Unit  Output  Unit 

[imperial  Unii  | 


□ 


r  ■; 

r 


Output  Options 

^  No  sreLvs_passes 

Material  Data  Output 

^  sref_vs_passes 

|No  Material  Output 

Analysis  Output  Options 
[Basic 

Number  of  Pass  to  Analyse 


Figure29:  The  Case  Control  tab 


Now  select  "File"  on  the  top  menu  and  select  "Save  As".  Save  the  input  file  in  the  same 
I  ocation  as  the  spectrum  fi  I  e  and  give  it  the  same  name  as  the  spectrum  fi  le.  N  ame  the  fi  I  e 
"Gen  seq .input",  see  Figure  30. 
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Figure  30:  T he  save  as  dialog  box 

The  input  f  i  I  e  I  ooks  si  mi  I  ar  to  that  presented  i  n  TabI  e  18. 
Table  18:  The  in  put  file 

Genseq -input 


Crack  Growth  Analysis  Program 
1 

-11 

1 

0 

1 

1 

7050-T7351 

2 

1 

0 

2 

2.5  5 

le-007 
2 

Flight  Hours:  100. 


7.3  Execution 

Oncethefiles  have  been  created,  click  on  "Build"  to  build  the  case,  see  Figure  31. 


File  Edit  View  Commands  Tools  Help 

i  -  H  X 

New  ;  Save  Remove 

Import  Export  Save  Out 

M  |■■■■■•  ♦ 

Build  Run 

\ 

Configure 

FigureSl:  CGAP  menu  bar 


You  should  see  a  new  label  underneath  cases  on  the  I  eft  hand  of  the  screen,  see  Figure  32. 
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File  Edit  View  Commands  Tools  Help 

l&  H  X 

New  Open  Save  Remove 

Ht 

Save  Out 

Build  Run 

ben_5eq.i:  Crack  Growth  Analysis 


V 


Material 


Load  I  Case  Cor 


Material  Database 

|7075TG51 

Material  Name  |DEC  1979  SAR-73-4 
Static  Properties 
Young's  Modulus  1 10300 


ksi 


Figure 32:  A  new  case  underneath  the  cases  control  box 


Multiplecasescan  be  loaded  into  CGAP  if  need  be. 


After  the  problem  has  been  successfully  built,  the  final  step  is  to  execute  the  analysis  by 
clicking  on  the  "Run"  button.  After  the  analysis  has  been  completed,  the  foil  owing  message 
should  bedisplayed  in  themessagewindow,  seeFigure33. 


Messages  |  Input  File  |  Output  File  |  Plot  File 


Spectrum  File  | 


Problem  loaded  successfullv.  Click  Run  to  analyse  it. 
Casel  running.. 

Analysis  Gen_seg.  input  completed  successfully. 


Figure 33:  M  essages  returned  after  analysis 


7.4  Results 

After  the  analysis  has  been  compi  eted  successful  ly  several  output  fi  I  es  are  created.  These  fi  I  es 
are  located  in  the  same  directory  as  the  input  file.  The  files  created  are: 

Gen_seq.output 

Gen_seq.passes 

Gen_seq.dmo 

Gen_seq.dsout 

Gen_seq.truncspc 

Gen_seq.passes  contains  the  notch  concentration  factors  and  their  associated  fatiguelifeto 
fai  lure  measured  i  n  fi  i  ght  hours.  For  this  exampi  ewe  expect  that  the  output  will  return  the 
val  ues  presented  i  n  TabI  e  19. 

Tablel9:  Theestimated  fatiguelifeat  each  notch  featureusing  CGAP 


Notch  concentration  factor 

Estimated  fatigue  iife(fiight  hours) 

2.5 

1,040,481 

5 

37,952 
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8.  Discussion 

Through  the  devd  opment  and  veri  fi  cati  on  of  the  CG  A  P  FA  M  5H  mod  ul  e  a  seri  es  of  test  cases 
were  developed.  These  cases  were  designed  to  test  certain  aspects  of  the  CGAP  FAMSFI 
module.  Someof  these  test  cases  have  been  included  in  Appendix  E  for  the  users  benefit.  A 
descri  pti  on  of  each  case  i  s  provi  ded  i  n  TabI  e  28.  These  test  cases  were  pri  mari  I  y  used  to  veri  fy 
the  correct  execution  of  the  program.  Itshould  denoted  thatthey  are  not  a  validation  of  the 
capability  of  FAMSFI  in  predicting  fatigueinitiation  lives.  All  test  cases  produced  consistent 
results  with  those  produced  by  the  FAM  SFI  executable. 

It  is  important  to  consider  the  limitations  with  the  current  strain-life  methodology.  An 
understanding  of  the  deficiencies  and  limitations  will  help  the  engineer  to  better  interpret 
resultsin  amoremeaningful  manner.  Following  isadiscussion  of  some  of  these  deficiencies 
and  limitations. 

8.1  Cyclic  Stress-Strain  Behaviour 

M  ateri  al  s  under  cycl  i  c  I  oad  i  ng  can  d  i  spl  ay  vari  ous  phenomena  that  are  not  as  yet  model  I  ed  i  n 
thecurrent  strain-life  algorithm  included  in  CGAP/  FAM  SFI.  The  FAM  SFI  module  evaluates 
damage  based  on  stabilised  stress-strain  hysteresis  loops.  But  it  is  well  known  that  this 
stabilised  stress-strain  behaviour  observed  undercyclicloading  may  bequitedifferentto that 
observed  under  monotonic  loading.  Obviously  this  has  significant  implications  in  the 
application  of  the  strain-life  approach.  For  accurate  estimates  of  strain  it  is  important  to  use 
the  correct  stress-strain  relationship.  The  difference  in  stress-strain  behaviour  between 
monotonic  loading  and  cyclic  loading  clearly  demonstrates  a  dependence  on  load  history. 
Thismaterial  dependence  on  thepriorloadingwas  first  observed  by  Bauschinger[29]and  is 
known  as  the  Bauschi  nger  effect.  Bauschinger  observed  that  the  yield  strength  in  tension  or 
compression  was  reduced  after  the  application  of  a  load  of  opposite  sign.  The  application  of 
one  single  load  is  enough  to  change  the  stress-strain  behaviour  of  a  material.  The  way  in 
which  the  behaviour  changes  after  one  single  load  is  also  affected  by  the  temper  of  the 
material.  A  dramatic  example  of  how  the  temper  of  a  material  can  alter  the  stress-strain 
relationship  with  repeated  loads  has  been  presented  by  Morrow  [30].  Fie  observed  the 
continuous  stress-strain  behaviour  of  copper  for  conditions  where  the  copper  was  (a)  fully 
annealed,  (b)  partially  annealed  and  (c)  cold  worked.  The  results  of  this  investigation  are 
presented  in  Figure  34.  For  this  particular  material  thecurves  eventually  stabilise  after  enough 
repeated  loads. 
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Figure  34:  Stress-strain  behavi  our  of  copper  subjected  to  cy  die  strain -controi  i  ed  axi  ai  ioads.  (a)Fuiiy 
anneaied,  showing  cydic hardening,  (b)  Partiaiiy  anneaied,  showing  smaii  cydic hardening 
and  softening,  (c)  Coid-worked,  showing  cydic  softening  [30]. 

Other  effects  such  as  mean  stress  relaxation  and  ratchetting  have  also  been  observed  with 
varying  degrees  in  different  materials.  Some  examples  of  the  how  these  effects  change  the 
stress-strain  responsewith  repeated  load  applications  are  provided  in  Figure  35.  Mean  stress 
relaxation  can  beobserved  under  strain-controlled  cyclic  loading  where  themean  stress  of  the 
spectrum  relaxes  tow  ardszero,  as  observed  in  Figure  35  (a).  When  anon-symmetric  constant 
amplitudeload  isapplied  under  stress  control  another  effect  known  as  ratchetting  may  occur 
as  in  Figure35(b).  At  structural  details  such  as  notches,  the  local  load  is  neither  under  strain 
control,  nor  under  stress  control,  irrespectiveof  the  mode  of  control  of  the  remote  loading. 
FI  ence,  we  may  expect  to  see  a  combi  nati  on  of  mean  stress  rel axati on  and  ratchetti  ng  as  i  n 
Figure  35(c). 


Figure35:  Eiastic-piastic  deformation  behaviour  subjected  to  (a)  constant  cydic  strains  (b)  constant 
cydic  stress^  and  (c)  remote  constant  nominai  cydic  stresses 
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In  FAMSH,  none  of  the  above  complex  stress-strain  behaviours  are  considered.  Instead,  the 
materi  al  s  are  assumed  to  be  of  M  asi  ng-type,  i.e. ,  the  hysteresis  I  oops  can  be  reconstructed  by 
magnifyi  ng  the  stead  y-statecycl  i  c  stress-strai  n  curve  by  a  factor  of  two.  Thesteady-statecycl  i  c 
stress-strain  curve  is  obtained  by  conducting  tests  on  smooth  (un-notched)  specimens 
subjected  to  fully-reversed  constant  amplitude  loading  [31]. 

Other  effects,  such  as  cyclic  softening  and  hardening,  ratchetting  etcwhich  effect  theshape  of 
the  hysteresis  loops,  have  not  as  yet  been  considered  i  n  the  code.  These  effects  are  hard  too 
quantify  in  situations  of  variable  amplitude  loading  and  require  a  more  sophisticated 
approach  when  modelling.  The  assumption  used  when  applying  the  FAMSH  code  is  that 
these  second  order  effects  have  minimal  impact  on  thefinal  solution.  However,  it  is  up  to  the 
user  to  ensure  the  validity  of  these  assumptions. 

8.2  N  otch  Root  Stress 

FA  M  SH  uses  the  N  euber  method  to  esti  mate  the  stress  at  the  root  of  a  notch.  The  N  euber 
method ,  w  h  i  I  e  assu  med  to  be  appi  i  cabi  e  to  a  w  i  d  e  range  of  geometri  es,  d  oes  have  I  i  mi  tati  ons. 
The  method  does  not  hold  true  in  situations  where  the  notch  root  stress  is  not  in  phase  with 
the  remote  I  oad .  It  does  not  account  for  ti  medependant  phenomenon  such  as  creep  and  stress 
relaxation.  It  also  does  not  account  for  cyclic  stress  relaxation  at  the  notch  root.  Whilethe 
Neuber  method  has  limitations  it  is  still  widely  used.  In  most  cases  the  Neuber  method 
sl  i  ghti  y  overesti  mates  the  notch-ti  p  stresses  and  strai  ns  [32].  Thi  s  may  resu  1 1  i  n  a  conservati  ve 
esti  mate  of  I  i  fe.  Other  methods  of  esti  mati  ng  the  notch  root  stress  are  avai  I  abl  esuch  as  that  by 
Glinka  [33],  which  has  been  shown  to  provide  better  results  in  conditions  of  planestrain  and 
torsional  loading  [34,  35].  The  inclusion  of  other  methods  into  CGAP  to  esti  mate  the  notch 
root  stress  could  improve  the  accuracy  of  predictions  for  certain  cases  where  the  Neuber 
method  has  been  shown  to  be  deficient.  At  present  the  N  euber  method  is  the  only  method 
available  and  it  must  be  remembered  that  the  method  may  lead  to  erroneous  solutions  in 
specific  cases. 

8.3  Definition  and  Validity  of  Material  Data 

Another  important  consideration  when  using  theFAMSH  codeisthe  origin  of  the  material 
data  used  i  n  the  analysis.  Thi  scan  be  demonstrated  by  consideri  ng  thei  nput  used  to  describe 
thestrain-lifecurve.Thiscurvecan  begenerated  in  several  ways  and  itisimportanttonoteon 
what  basisthese  curves  were  created.  A  key  consideration  isthedefinition  of  failure.  Failure 
could  be  considered  to  haveoccurred  when  the  test  specimen  reaches  a  predefined  level  of 
compliance  or  when  thespeci  men  has  completely  failed.  A  changeinthefai  I  uredefiniti  on  can 
completely  change  the  interpretation  of  the  solution.  Another  consideration  is  the  surface 
conditi  on  of  thetest  speci  mens  used  to  generate  thestrai  n-l  ifecurve.  T  ypi  cal  I  y  these  tests  are 
performed  with  smooth  (un-notched)  specimens,  with  certain  amount  of  surface  polishing, 
but  the  su  rface  cond  i  ti  on  of  the  structu  re  may  not  al  w  ays  match  those  of  the  speci  mens.  A  gai  n 
this  will  affect  interpretation  of  the  results  produced  by  FAMSH  and  the  engineer  must 
consi  d  er  th  i  s  w  hen  i  nterpreti  ng  the  resu  I  ts. 

It  is  worth  highlightingthestrain-life  techniques  sensitivity  to  small  strai  ns  wherefatiguelife 
estimates  are  high.  To  demonstrate  this  issue,  consider  two  models  both  using  the  same 
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material,  oneusingtheequati  on  description  and  the  other  usingthetabu  I  ar  input.  The  models 
areidentical  and  onlydiffer  inthewaythenotionally-identical  material  properties  are  entered 
into  the  program.  After  running  the  same  problem  using  both  the  material  inputs,  the 
following  output  is  produced,  see  Table  20. 

Table20:  Total  life  predicted  for  various  K„  values 


Kn 

Equation 

Input 

Tabular 

Input 

%  Difference 

3 

153341.8 

155716.2 

-1.52% 

3.5 

77144.6 

78133.6 

-1.27% 

4 

40037.6 

40235.6 

-0.49% 

4.5 

21161.3 

21903.5 

-3.39% 

5 

13673.3 

14215.3 

-3.81% 

5.5 

9718.0 

9004.4 

7.92% 

6 

5632.5 

5748.5 

-2.02% 

6.5 

3786.7 

3693.3 

2.53% 

7 

2449.8 

2321.4 

5.53% 

Although  the  inputs  for  the  tabular  definition  match  the  equation  definition  exactly,  the 
results  in  Table  20  indicatethat  errors  introduced  during  interpolation  have  contributed  to 
d i screpanci  es  i  n  thef i  nal  resu I ts.  The associ  ated  sensi ti vi ty  to  i  n put  d ata  parti  cu I  arl y  for  smal  I 
strains  and  long  lives  should  betaken  intoconsideration  when  interpreting  strain-liferesults 
using  thistechnique. 

8.4  Damage  Summation 

TheFAMSH  code  uses  Miner's  [25]  ruletodeterminethe  level  of  damage.  Miner's  rule  uses 
the  linear  accumulation  concept  proposed  by  Palmgren  [26]  as  a  measure  of  damage.  This 
method  of  evaluating  the  amount  of  damaged  incurred  by  each  cycle  does  not  account  for 
load  hi  story  effects.  In  other  words.  Miner's  rule  does  not  consider  the  effect  the  sequence  of 
I  oads  has  on  the  accumul  ati  on  of  damage  i  n  a  structure.  I  gnori  ng  these  effects  can  result  i  n  an 
order  of  magnitude  difference  in  thefi  nal  result  [36, 37].  Improvements  to  the  code  could  be 
implemented  via  the  introduction  of  more  sophisticated  non-linear  damage  accumulation 
algorithms.  A  number  of  such  models  can  be  found  in  the  literature  [27]. 
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9.  Future  Improvements 

The  fol  I  ow i  ng  mod i f i  cati  ons  to  the  FA  M  SH  cod e  are  u nder  consi  d erati  on  for  f utu re  C G A  P 
releases. 

•  The  incorporation  of  a  cyclic  strain  hardening  model 

•  Ability  to  queue  and  run  multiplejobs  (similar  to  a  batch  file  process) 

•  Additional  and  improved  equivalent  strain  equations 

•  Moreadvanced  damage  accumulation  modelsto  account  for  sequence  effects 

•  The  i  ncl  usi  on  of  other  techni  ques  to  esti  mate  the  notch  root  stress 

•  Seamless  interaction  between  strain-life  and  crack  growth  algorithms  to  providetotal 
life  esti  mates. 

•  Inclusion  or  emulation  of  other  strain-life  codes 

Medium  to  long  term  research  activities  that  may  benefit  from  the  CGAP  analysis 
envi  ronment  and  i  mprove  the  anal  ysis  capabi  I  ity  i  ncl  ude  a  study  of  the  effect  of  net  secti  on 
size  and  surface  condition  on  the  fatigue  limit.  Improved  damage  models  may  help  to 
disassociate  the  portion  of  damage  due  to  crack  growth  and  the  component  due  to  strain 
harden!  ng  of  the  materi  al .  Thi  s  may  I  ead  to  more  robust  materi  al  properti  es  that  may  I  ead  to  a 
generalised  predictivecapability.  More  recent  developments  in  fatigue  life  predict! on  have 
utilised  the  strain-life  approach  to  predict  crack  growth  directly  [38-41].  Further  research  and 
development  of  this  approach  could  be  introduced  to  the  CGAP  environment  and  allow 
further  evaluation  of  this  techni  qua 


10.  Conclusion 

Thi  s  docu  ment  su  mmari  ses  the  i  mportant  aspects  of  the  strai  n-l  ife  theory  and  expl  ai  ns  how 
th  i  s  method  i  s  i  mpl  emented  i  nto  the  CG  A  P  envi  ronment  usi  ng  the  strai  n-l  i  fe  mod  u  I  e  FA  M  SFI . 
The  document  also  provides  the  reader  with  guidance  on  and  presents  examples 
demonstrating  the  use  of  the  new  module. 

The  ad  d  i  ti  on  of  the  FA  M  SFI  strai  n- 1  i  f  e  al  gori  th  m  i  nto  the  CG  A  P  tool  set  i  ntrod  uces  consi  stent 
materi  al  data  hand  I  i  ng,  bui  I  t-i  n  error  checki  ng,  and  a  fami  I  i  ar  graphi  cal  user  i  nterface,  ai  d  i  ng 
quality  management  processes.  Theseimprovementswill  I  ead  to  improved  robustness  in  the 
pred  i  cti  ons  w  hi  ch  ul  ti  mate!  y  contri  butes  to  more  effi  ci  ent  and  rel  i  abl  e  management  of  ai  rcraft 
structural  integrity. 
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Appendix  A;  The  Basis  of  the  FA  MSH  Stress-Strain 

Reiationship 

In  FAMSH,  the  following  equation  is  used  to  relate  the  plastic  stress  and  the  plastic  strain 
amplitude, 

^  >  (^pl  (27) 

where  c7^  is  the  plastic  stress,  is  the  proportional  limit,  isthetotal  applied  stress  and 

K'  and  n'  are  the  cyclic  strength  coefficient  and  cyclic  strain  hardening  coefficients 
respectively. 

U  si  ng  si  mi  I  ar  assumpti  ons  to  that  of  the  Ramberg-Osgood  rel  ati  on,  we  can  express  the  total 
strain  as  the  sum  of  the  elastic  and  plastic  components  of  strain  such  that: 


Sa^£e+^p 

(28) 

Rearranging  equation  27,  we  get 

/ 

1  / 

V 

(29) 

Therefore  equation  28  becomes 

(30) 

^  (^pl 

F  —  — — 

E 

Equation  30  can  be  used  to  represent  the  stabilised  hysteresis  I  oops  of  the  stress-strain  relation 
ship  using  a  Masing  type  material  approximation. 
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AppendixB;  Deriving  Material  Parameters  for  the 
FAMSH  Strain  Life  Equation 

The  strain-life  equation  used  in  FAMSH  hasahighdegreeof  flexibility,  however  determining 
material  parameters  for  this  equation  is  not  immediately  apparent.  This  section  provides  a 
brief  overview  of  one  technique  that  can  be  used  to  determine  the  input  parameters  for  the 
FAMSH  strain-life  equation.  The  method  presented  here  is  one  method,  but  the  reader  is 
remi  nded  that  there  maybe  other  better  methods  of  eval  uati  ng  these  parameters. 

Forthepurposesof  thisguide,  wehavechosen  to  usethetabulated  data  for  7050-T7351from 
the CG A P  material  database.  By  using  thetermsPLET  andELET  wecandividethestrain-life 
curveintothreesections,seeFigure36.TheparametersPLET  and  ELET  aredetermined  bythe 
user  and  the  avail  able  tabulated  data. 


Region  3 


Using  the  tabulated  data  for  the  material  7075-T7351  wecan  create  three  data  sets.  The  first 
data  set  contains  all  thedatawith  astrain  abovePLET  (thatisall  thedatain  Region  l).The 
second  data  set  contains  thedata  in  Region  2  between  PLET  and  ELET,  and  thethird  dataset 
contains  all  the  data  with  a  strain  below  ELET  (Region  3).  After  separating  the  7075-T7351 
data  into  these  three  sets,  the  tabulated  data  should  look  something  like  that  in  Table  21 
depending  on  your  choice  of  PLET  and  ELET. 
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Table21:  7075-T7351  material  data  separated  into  threesets 


Original  Data  Set 

Data  Set  1 
[Region  1) 

Data  Set  2 
[Region  2) 

DataSets 
[Region  3) 

Life 

Strain 

Life 

Strain 

Life 

Strain 

Life 

Strain 

10 

0.085 

10 

0.085 

20 

0.051 

20 

0.051 

30 

0.039 

30 

0.039 

40 

0.0325 

40 

0.0325 

60 

0.025 

60 

0.025 

88 

0.02 

88 

0.02 

145 

0.0155 

145 

0.0155 

200 

0.013 

200 

0.013 

350 

0.0105 

350 

0.0105 

400 

0.01 

400 

0.01 

500 

0.0093 

500 

0.0093 

600 

0.0087 

600 

0.0087 

700 

0.0083 

700 

0.0083 

900 

0.0078 

900 

0.0078 

1500 

0.007 

1500 

0.007 

2000 

0.0067 

2000 

0.0067 

3000 

0.0063 

3000 

0.0063 

5000 

0.0056 

5000 

0.0056 

40000 

0.0032 

40000 

0.0032 

1000000 

0.0015 

1000000 

0.0015 

In  this  example  PL  ET  =0.0105and  ELET  =0.0063.Thereare8other  unknown  parameters  to 
be  determined,  and  they  areCl,  C2,  Bl,  Gl,  C3,  C4,  B2and  G2.  Perhaps  the  easiest  method  to 
determi  nethese  parameters  is  to  use  a  non-l  i  near  curvefitti  ng  program  such  as  that  contai  ned 
in  Origin® 6.0  [42], 

If  you  use  Origin®  6.0,  you  will  need  to  create  the  function  used  in  the  fitting  process.  The 
inputs  should  look  something  like  that  shown  in  Figure  37. 
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Figures?:  U  ser  defined  non-iinear  curve  fitting  function 

U  si  ng  the  non-l  i  near  curvefi  tti  ng  program  determi  nethe  parameters  C 1,  C4,  B1  and  G2  usi  ng 
the  second  data  set.  As  we  have  chosen  PLET  to  be  greater  than  ELET,  Region  2  is  defined 
with  the  parameters  Cl,  C4,  Bland  G2  (If  PLET  was  less  than  ELET  then  Region  2  would  be 
defined  with  C2,C3,  Gland  B2).  The  fit  should  look  something  likethatshown  in  Figure  38. 


Figure 38:  C urve  fit  for  data  set  2  using  0  rigin®  6.0 

When  interpreting  the  results  of  the  non-linear  curvefit  it  is  necessary  to  identify  the  terms 
relating  to  the  elastic  strain  and  the  plastic  strain.  Remember  that  the  strain-life  equation 
(equation  12)  isthesumoftheplasticand  elastic  components  of  strain.  The  plastic  component 
of  the  equati  on  i  s  i  dentifi  ed  by  the  I  arger  grad  i  ent  on  a  I  og-l  og  pi  ot  Thus  the  I  arger  of  the  two 
exponents  in  equation  12  identifies  the  plastic  component. 
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By  fitting  Region  2  first  we  now  have  values  for  Cl,  C4,  B1  and  G2  which  we  can  use  in 
determiningtheparametersintheother  regions.  Now  usingthefirst  data  set,  inputCland  B1 
as  fixed  parameters  and  use  the  non-linear  curve  fitting  algorithm  to  find  C2  and  Gl.  The 
result  of  thefitting  should  look  likethat  shown  in  Figure  39. 


Similarly  we  now  use  C4  and  G2  (determined  from  our  non-linear  curve  fit  of  Region2)  as 
fixed  parameters  to  determi  nethefi  nal  two  unknown  parameters  C3  and  B2  usi  ng  data  set  3. 
Thecurvefit  is  presented  in  Figure 40. 


Figure  40:  C urve  fit  for  data  set  3  using  0  rigin  ®  6.0 

After  fittingthefi  nal  setof  data,wenow  have  deter  mined  all  the  parameters  needed  to  define 
the  strai  n-l  ife  equati on.  These  parameters  are  presented  i  n  TabI  e  22. 
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Table22:  Fitting  parameters  for  theFA  M  SH  strain-iife equation  for  7075-T7351  aiuminium 


Region  2 


Region  1 

Region  3 

Parameter 

C2 

G1 

Cl 

B1 

C4 

G2 

C3 

B2 

Value 

0.5087 

0.81999 

0.00907 

0.05491 

1.37415 

0.99716 

0.04192 

0.24271 

In  this  example  when  compared  to  the  original  7075-T7351  aluminium  data,  the  FAMSH 
strain-life  equation  produces  excel  lent  agreement,  see  Figure  41. 


Figure41:  TheFAMSH  strain-iife  equation  compared  to  the  origin  ai  7075-T7351  aiuminium  data 
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AppendixC;  FAMSH  OutputFiles 

Table23:  An  exampleofthefilenamedmofile 


Example.dmo 

FAMSH  PROGRAM  INITIATED  20090731  10:45:06 

Title:  generic  test  seq 

LOADS  APPLICATION  NUMBER  TWO 

KN  =  1.0000 

SREF  =  2.500 

(KN*SREF)  = 

2.500 

DAMAGE  DISTRIBUTION 

BY  MISSION 

MISSION 

MISSION 

TOTAL 

AVG  DMG 

%  MISSION 

NUMBER 

UTILIZ . 

DAMAGE 

PER  FLT .  DMG 

OF  TOTAL 

1 

121 

0 . 1345E-02 

0 . llllE-04 

0.88 

2 

65 

0 . 1545E-02 

0.2377E-04 

1.01 

3 

63 

0.2486E-02 

0.3946E-04 

1.63 

4 

132 

0 . 1991E-02 

0.1508E-04 

1.30 

5 

94 

0.2596E-02 

0.2762E-04 

1.70 

6 

139 

0 . 6220E-02 

0.4475E-04 

4.07 

7 

75 

0 . 1867E-02 

0.2489E-04 

1.22 

8 

193 

0.5910E-02 

0.3062E-04 

3.87 

9 

96 

0.3295E-02 

0.3433E-04 

2.16 

10 

95 

0.2777E-02 

0.2923E-04 

1 . 82 

11 

114 

0.4070E-02 

0.3570E-04 

2.66 

12 

33 

0.6235E-03 

0 . 1889E-04 

0.41 

13 

91 

0.3849E-02 

0.4230E-04 

2.52 

14 

126 

0.4775E-02 

0.3789E-04 

3.12 

15 

111 

0.2781E-02 

0.2505E-04 

1 . 82 

16 

175 

0.6870E-02 

0.3926E-04 

4.49 

17 

216 

0.1363E-01 

0.6309E-04 

8.91 

18 

206 

0 . 1511E-01 

0.7336E-04 

9.88 

19 

77 

0 . 9120E-02 

0 . 1184E-03 

5.96 

20 

102 

0.4437E-02 

0.4350E-04 

2.90 

21 

66 

0.4884E-02 

0.7399E-04 

3.19 

22 

82 

0.5841E-02 

0.7124E-04 

3.82 

23 

153 

0 . 1823E-02 

0 . 1192E-04 

1.19 

24 

107 

0.2016E-02 

0 . 1884E-04 

1.32 

25 

87 

0 . 1791E-02 

0.2059E-04 

1 . 17 

26 

90 

0.2324E-02 

0.2582E-04 

1.52 

27 

120 

0.3334E-02 

0.2778E-04 

2 .18 

28 

109 

0.4402E-02 

0.4038E-04 

2.88 

29 

102 

0.4158E-02 

0.4077E-04 

2.72 

30 

191 

0.3698E-02 

0.1936E-04 

2.42 

31 

143 

0.1063E-02 

0.7436E-05 

0.70 

32 

94 

0 . 1027E-02 

0.1093E-04 

0.67 

33 

290 

0.2573E-02 

0.8874E-05 

1.68 

34 

150 

0.5210E-02 

0.3473E-04 

3.41 

35 

127 

0 . 1492E-02 

0 . 1175E-04 

0.98 

36 

100 

0 . 1847E-02 

0 . 1847E-04 

1.21 

37 

55 

0 . 9397E-02 

0.1709E-03 

6.15 

38 

11 

0.7360E-03 

0 . 6691E-04 

0.48 

TOTAL 

4401 

0.1529 

0.3475E-04 

100.0 
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T  able24:  A  n  example  of  the  damage  density  information  contained  in  the  filenamedsout  file 


Example.dsout 

FAMSH 

PROGRAM 

INITIATED 

20090604 

Title: 

generic  test  seq 

LOADS 

APPLICATION  NUMBER 

TWO 

KN  = 

1.0000 

SREF  = 

2.500 

(KN*SREF)  = 

2.500 

OCCURRENCES 

OCCURRENCES 

BIN 

#  RANGE  (%)  IN 

RANGE 

PCT  DMG 

REMAINING 

1 

0  - 

5 

14.5 

0.0 

441.5 

2 

5  - 

10 

54.0 

1.7 

387.5 

3 

10  - 

15 

347.0 

41 . 6 

40.5 

4 

15  - 

20 

34.0 

22 . 8 

6.5 

5 

20  - 

25 

4.5 

5.6 

2.0 

6 

25  - 

30 

0.5 

1.8 

1.5 

7 

30  - 

35 

0.0 

0.0 

1.5 

8 

35  - 

40 

0.0 

0.0 

1.5 

9 

40  - 

45 

0.5 

7.3 

1.0 

10 

45  - 

50 

1.0 

19.1 

0.0 

11 

50  - 

55 

0.0 

0.0 

0.0 

12 

55  - 

60 

0.0 

0.0 

0.0 

13 

60  - 

65 

0.0 

0.0 

0.0 

14 

65  - 

70 

0.0 

0.0 

0.0 

15 

70  - 

75 

0.0 

0.0 

0.0 

16 

75  - 

80 

0.0 

0.0 

0.0 

17 

80  - 

85 

0.0 

0.0 

0.0 

18 

85  - 

90 

0.0 

0.0 

0.0 

19 

90  - 

95 

0.0 

0.0 

0.0 

20 

95  - 

100 

0.0 

0.0 

0.0 

21 

100 

0.0 

0.0 

0.0 

TOTAL  CYCLES 

456.0 

TOTAL  DAMAGE 

1.22661E-04 

1  MAXIMUM  DELTA  STRESS: 

60000.00 

T able  25:  A  n  example  of  the  filenamepasses  file 


Example.passes 

C:\Cgap\Example  l\Example. 

Flight  Hours 

Kn 

8152,52148437500 

2.500000 

313,114379882812 

5.000000 

Table26:  An  exampleof  the filenametruncspc  file 


Example.truncspc 

2 

1 

CYCLE 

0 

1.000 

CYCLE  CCUNT  &  FLIGHT 

DURATICN 

4391 

38 

FLIGHT  NUMBER  &  MISSICN  TYPE 

224 

1.000 

CYCLE  CCUNT  &  FLIGHT 

DURATION 

0 

0.00 

38001022 

1.00 

38001022 

-13408.00 

38002023 

-1773.00 

38002023 

-9068.00 

38003019 

6891.00 

38003019 

1087.00 

38004004 

15104.00 

38004004 

5740.00 

38004005 

13942.00 

38004005 

10422.00 

38004004 

13439.00 

38004004 

7405.00 

38005004 

13167.00 

38005004 

7246.00 

38007008 

13512.00 

38007008 

10301.00 

38007007 

14201.00 

38007007 

6401.00 

38007007 

14201.00 

38007007 

6401.00 

38007007 

13380.00 

38007007 

7222.00 

38007007 

13380.00 

38007007 

7222.00 

38007007 

13380.00 

38007007 

7222.00 

38007007 

14714.00 
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Table27:  An  exampleof  the  filenameoutput  file 


Example-output 


***  FAMSH  Version  1.42  -  Modified  FAMS 

FAMS  PROGRAM  INITIATED  20090714  13:24:17 

*  *  *  *  Information  *  *  *  * 

Passes  vs  Kn  data  are  appended  to  file 

C:\Program  f iles\DSTO\CGAP\Analysis  Files\Example . passes 


*****************  Information  *****-* 
Default  (Morrow)  equation  selected. 


FATIGUE  ANALYSIS  OE  METALLIC  STRUCTURES 
(FAMS) 


Developed  by  N.R.Krishnan 
Modified  by  K.N. Bailey 


MATERIAL  DATA  Material 

Young's  Modulus  (psi) 
Yield  Stress  (psi) 
Ultimate  Stress  (psi) 
Fracture  Stress  (psi) 


D7075-T651  AL  ALLOY  (DSTO  Trendline) 
1.0300E+07 
68000. 

82500. 

110000. 


Title:  generic  test  seq 

LOADS  APPLICATION  NUMBER  TWO 

KN  =  1.0000  SREF  =  2.500  {KN*SREF)  =  2.500 

name  2 

TIME  REPRESENTED  BY  LOADING  FOR  BLOCK  NO .  1  IS  1.000 

NOTCH  STRESSES  AND  STRAINS  COMPUTED  FOR  458  UNIQUE  LEVELS 
BLOCK  DAMAGE  =  6.78602E-05 


name  2 

TIME  REPRESENTED  BY  LOADING  FOR  BLOCK  NO .  2  IS  1.000 

NOTCH  STRESSES  AND  STRAINS  COMPUTED  FOR  456  UNIQUE  LEVELS 
BLOCK  DAMAGE  =  5.48013E-05 

TOTAL  DAMAGE  =  1.22661E-04  LIFE  =  16305.0  NUMBER  PASSES  =  8152.5215 

APPLIED  MAXIMUM  LOAD  =  16952.00  MINIMUM  LOAD  =-13408.00 

NUMBER  OE  BLOCKS  <ELIGHTS>  ANALYZED  =  2 

LOADS  APPLICATION  NUMBER  TWO 

name  2 

NOTCH  STRESSES  AND  STRAINS  COMPUTED  FOR  458  UNIQUE  LEVELS 
BLOCK  DAMAGE  =  6.78602E-05 

name  2 


NOTCH 

STRESSES  AND 

STRAINS  COMPUTED  FOR 

456  UNIQUE  LEVELS 

BLOCK 

DAMAGE  =  5 . 

48013E-05 

TOTAL  DAMAGE  = 

1 

22661E-04 

LIFE  = 

16305.0 

NUMBER  PASSES 

=  8152. 

5215 

LI EE  UNDER  LOAD 

BLOCKS  REPEATED 

16305.0 

NUMBER  PASSES 

=  8152. 

5215 

The  option  was 

chosen  to  use 

LOADS  APPLICATION  NUMBER  TWO  for  damage 

and  truncation 

studies . 

TRUNCATION 

%  OE  ORIGINAL 

CYCLES 

CYCLES 

% 

CUTOEF 

DAMAGE 

DAMAGE 

REMOVED 

REMAINING 

REMOVED 

l.OOE-08 

0.00012 

100.00 

30.00 

426.00 

6 

.579 

l.OOE-07 

0.00011 

89.38 

192.00 

264.00 

42 

.105 

l.OOE-06 

0.00005 

39.07 

442.00 

14.00 

96 

.930 

l.OOE-05 

0.00004 

28.54 

451.00 

5.00 

98 

.904 

l.OOE-04 

0.00004 

28.54 

451.00 

5.00 

98 

.904 

l.OOE-03 

0.00004 

28.54 

451.00 

5.00 

98 

.904 

l.OOE-02 

0.00004 

28.54 

451.00 

5.00 

98 

.904 

l.OOE-01 

0.00004 

28.54 

451.00 

5.00 

98 

.904 

l.OOE-07 

0.00011 

89.38 

192.00 

264.00 

42 

.105 

TOTAL  #  CYCLES 

456.00 

TOTAL  #  REAL 

CYCLES 

455.00 

TOTAL  #  VIRTUAL 

CYCLES  'RAINELOW': 

1.00 

TOTAL  #  CYCLES 

456.00 
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RESIDUAL  STRESS  AT  END  OF  LOADING  =  0.00 

LAST  LOAD  CAUSING  YIELD  =  0.00 


Title:  generic  test  seq 

LOADS  APPLICATION  NUMBER  TWO 

KN  =  1.0000  SREF  =  5.000  {KN*SREF)  =  5.000 

name  2 

TIME  REPRESENTED  BY  LOADING  FOR  BLOCK  NO .  1  IS  1.000 

NOTCH  STRESSES  AND  STRAINS  COMPUTED  FOR  458  UNIQUE  LEVELS 
BLOCK  DAMAGE  =  1.87001E-03 

name  2 

TIME  REPRESENTED  BY  LOADING  FOR  BLOCK  NO .  2  IS  1.000 

NOTCH  STRESSES  AND  STRAINS  COMPUTED  FOR  456  UNIQUE  LEVELS 
BLOCK  DAMAGE  =  1.51884E-03 

TOTAL  DAMAGE  =  3.38886E-03  LIFE  =  590.169  NUMBER  PASSES  =  295.0847 

APPLIED  MAXIMUM  LOAD  =  16952.00  MINIMUM  LOAD  =-13408.00 

NUMBER  OF  BLOCKS  <FLIGHTS>  ANALYZED  =  2 


LOADS  APPLICATION  NUMBER  TWO 


name  2 

NOTCH  STRESSES  AND  STRAINS  COMPUTED  FOR  458  UNIQUE  LEVELS 
BLOCK  DAMAGE  =  1.68650E-03 


name  2 

NOTCH  STRESSES  AND  STRAINS  COMPUTED  FOR  456  UNIQUE  LEVELS 
BLOCK  DAMAGE  =  1.50659E-03 

TOTAL  DAMAGE  =  3.19310E-03  LIFE  =  626.351  NUMBER  PASSES  =  313.1757 

LIFE  UNDER  LOAD  BLOCKS  REPEATED  =  626.229  NUMBER  PASSES  =  313.1144 

The  option  was  chosen  to  use  LOADS  APPLICATION  NUMBER  TWO  for  damage  and  truncation  studies. 


TRUNCATION 

%  OF  ORIGINAL 

CUTOFF 

DAMAGE 

DAMAGE 

l.OOE-08 

0.00319 

100.00 

l.OOE-07 

0.00319 

100.00 

l.OOE-06 

0.00318 

99.68 

l.OOE-05 

0.00169 

52.91 

l.OOE-04 

0.00039 

12.11 

l.OOE-03 

0.00039 

12.11 

l.OOE-02 

0.00039 

12.11 

l.OOE-01 

0.00039 

12.11 

l.OOE-07 

0.00319 

100.00 

TOTAL  #  CYCLES 

TOTAL  # 

REAL  CYCLES 

:  455.00 

TOTAL  # 

VIRTUAL 

CYCLES 

RAINFLOW':  1.00 

TOTAL  # 

CYCLES 

:  456.00 

RESIDUAL 

STRESS  AT 

END  OF 

LOADING  =  -3561 . 05 

LAST  LOAD 

CAUSING 

YIELD 

16952.00 

CYCLES 

CYCLES 

% 

REMOVED 

REMAINING 

REMOVED 

10.00 

446.00 

2.193 

10.00 

446.00 

2.193 

29.00 

427.00 

6.360 

394.00 

62.00 

86.404 

451.00 

5.00 

98.904 

451.00 

5.00 

98.904 

451.00 

5.00 

98.904 

451.00 

5.00 

98.904 

10.00 

456.00 

446.00 

2.193 

Analysis  completed  successfully. 
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Appendix  D;  Sample  FAMSH  Materials  Files 


Sample  f/matl.dat  file 

>D7075-T651  AL  ALLOY  (DSTO  Trendline) 

10300  68  82.5  110 

32  1000  1 

0.0058732  60.494 

0.0063586  65.493 

0.0068504  70.427 

0.0074646  74.137 

0.0082181  76.542 

0.0090608  78.314 

0.0099697  79.743 

0.010935  80.959 

0.0119517  82.031 

0.013017  82.999 

0.0141289  83.886 

0.0152862  84.711 

0.0164878  85.485 

0.0177327  86.216 

0.0190204  86.911 

0.0203503  87.576 

0.0217217  88.214 

0.0231342  88.828 

0.0245875  89.421 

0.0260812  89.996 

0.0276152  90.554 

0.0291896  91.096 

0.0308049  91.625 

0.0324628  92.142 

0.0341662  92.649 

0.0359213  93.149 

0.0377396  93.644 

0.0396416  94.141 

0.0416642  94.646 

0.0438716  95.176 

0.053135  97.179 

0.0654031  99.435 

20 

0.018  602 
0.013  1100 

0.012  1340 

0.0085  2980 

0.008  3640 

0.0075  4450 

0.007  6000 

0.0065  7330 

0.006  9900 

0.0055  12100 

0.005  14800 

0.0045  22000 

0.004  32900 

0.0035  54200 

0.003  109000 

0.0025  243000 

0.002  540000 


Sample  f/matZdat  file 

>7075-1651  DEC. 1979  SAR-79-4 

10300  68  82.5  110  1.31575e-005  3.27 

0.01103  0.3221  0.0035  2.56 

0.0588  0.1573  0.7444  3.1484  0.01151  0.009454 
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AppendixE:  FAMSH  Test  Cases 

T 0  veri  fy  the  FA  M  5H  mod  u  I  e,  a  n  u  mber  of  test  cases  w  ere  d  evel  oped .  Some  of  the  test  cases 
have  been  included  here.  All  the  test  cases  were  evaluated  using  the  original  FAMSFI  code 
and  compared  to  the  results  produced  byCGAP.AII  cases  returned  consistent  results.  Note: 
These  mod  el  s  are  on  I  y  to  be  used  for  val  I  d  ati  on  and  not  for  any  other  pu  rposes.  M  ateri  al  d  ata 
and  parameters  should  not  be  transferred  tosolvereal  engineering  problems. 

The  attached  CD  ROM  contains  all  theCGAP  strain  life  module  val  id  ati  on  cases. 

Table28:  FAMSH  Test  Casedescriptions 

Case  No.  Test  Description 

1  Uses  a  moderate  si  zed  spectrum  file  in  conjunction  with  two  reference  stress  values 

2  This  case,  similar  to  the  first,  uses  a  smaller  version  of  the  same  spectrum  file  this  time  only 
containing  2  flights 

3  This  case  uses  four  reference  stresses  and  isdesigned  to  fail  on  the  last  reference  stress  analysed 

4  A  spectrum  file  containing  2  flights  is  defined  in  block  format  using  DSA 

5  A  spectrum  file  containing  1  flight  is  defined  in  block  format  using  DSA 

6  A  spectrum  file  containing  1  flight  is  defined  in  cycle  format  using  DSA 

7  A  spectrum  file  containing  1  flight  is  defined  in  cycleformat  with  no  DSA 

8  A  spectrum  file  containing  2  flights  is  defined  in  cycleformat  with  no  DSA 

9  A  spectrum  file  containing  1  flight  is  defined  in  block  format  with  no  DSA 

10  A  spectrum  file  containing  2  flights  is  defined  in  block  format  with  no  DSA 

11  A  large  spectrum  file  is  used 

12  An  equation  material  description  is  used  (7075-T651)  with  the  case  1  spectrum 

13  Removed  the  Flight  hours  and  maximum  delta  stress  option  from  the  input  file 

14  No  stress-strain  information  is  output  and  Basic  analysis  output  is  requested 

15  Cyclic  stress-strain  information  is  output  and  Basic  analysis  output  is  requested 

hysteresis  stress  versus  hysteresis  stress-strain  information  isoutputand  Basic  analysis  output  is 
requested 

17  No  stress-strain  information  is  output  and  detailed  analysis  output  is  requested 

18  No  stress-strain  information  isoutputand  moredetailed  analysis  output  is  requested 

19  No  stress-strain  information  is  output  and  most  detailed  analysis  output  is  requested 

20  Residual  stress  is  added  to  the  analysis 

21  One  pass  is  used  in  the  analysis 

22  A  Kn  of  1.5is  used  in  conjunction  with  the  reference  stresses 

23  A  conversion  factor  of  1.2  is  used  in  the  analysis 

24  Theloopin  equivalent  strain  equation  is  requested 

25  The  modified  loopin  equivalent  strain  equation  is  requested 

26  The  Walker  equivalent  strain  equation  is  requested 

27  TheSmith-Watson-Topper  equivalent  strain  equation  is  requested 

28  The  F-18  equivalent  strain  equation  is  requested _ 
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Table29:  FAMSH  Test  Caseinputs 


Case  No. 

IPTl 

IPT2 

IPT3 

IPT4 

IPT5 

Kn 

Material 

MATF  UNITCV 

RS 

NMLT 

SREF 

DCUT 

ESE 

m  OC 

>0 

y  Flight 
’  Hours 

DSA 

CYCLE/'^“'’';P'^ 

—  %'S: 

1 

1 

12 

1 

0 

0 

1 

D7075-T651 

1 

1 

0 

2 

2.5 

5 

l.E-07 

2 

15000 

YES 

CYCLE  Multiple 

2 

1 

12 

1 

0 

0 

1 

D7075-T651 

1 

1 

0 

2 

2.5 

5 

l.E-07 

2 

15000 

YES 

CYCLE  Multiple 

3 

1 

12 

1 

0 

0 

1 

D7075-T651 

1 

1 

0 

4 

2.5 

5  1  10 

l.E-07 

2 

15000 

YES 

CYCLE  Multiple 

4 

1 

12 

1 

0 

0 

1 

D7075-T651 

1 

1 

0 

2 

2.5 

5 

l.E-07 

2 

15000 

YES 

BLOCK  Single 

5 

1 

12 

1 

0 

0 

1 

D7075-T651 

1 

1 

0 

2 

2.5 

5 

l.E-07 

2 

15000 

YES 

BLOCK  Multiple 

6 

1 

12 

1 

0 

0 

1 

D7075-T651 

1 

1 

0 

2 

2.5 

5 

l.E-07 

2 

15000 

YES 

CYCLE  Single 

7 

1 

12 

1 

1 

0 

1 

D7075-T651 

1 

1 

0 

2 

2.5 

5 

l.E-07 

2 

15000 

NO 

CYCLE  Single 

8 

1 

12 

1 

1 

0 

1 

D7075-T651 

1 

1 

0 

2 

2.5 

5 

l.E-07 

2 

15000 

NO 

CYCLE  Multiple 

9 

1 

12 

1 

1 

0 

1 

D7075-T651 

1 

1 

0 

2 

2.5 

5 

l.E-07 

2 

15000 

NO 

BLOCK  Single 

10 

1 

12 

1 

1 

0 

1 

D7075-T651 

1 

1 

0 

2 

2.5 

5 

l.E-07 

2 

15000 

NO 

BLOCK  Multiple 

11 

O 

-21 

1 

0 

0 

1 

D7075-T651 

1 

1 

0 

2 

2.5 

5 

l.E-07 

2 

15000 

YES 

CYCLE  Multiple 

12 

1 

2D 

1 

0 

0 

1 

7075-T651 

2 

1 

0 

2 

2.5 

5 

l.E-07 

2 

15000 

YES 

CYCLE  Multiple 

13 

1 

U. 

1 

0 

0 

1 

D7075-T651 

1 

1 

0 

2 

2.5 

5 

l.E-07 

2 

1 

YES 

CYCLE  Multiple 

14 

o 

-21 

1 

0 

0 

1 

D7075-T651 

1 

1 

0 

2 

2.5 

5 

l.E-07 

2 

15000 

YES 

CYCLE  Multiple 

15 

1 

-21 

1 

0 

0 

1 

D7075-T651 

1 

1 

0 

2 

2.5 

5 

l.E-07 

2 

15000 

YES 

CYCLE  Multiple 

16 

2 

-21 

1 

0 

0 

1 

D7075-T651 

1 

1 

0 

2 

2.5 

5 

l.E-07 

2 

15000 

YES 

CYCLE  Multiple 

17 

O 

2D 

1 

0 

0 

1 

D7075-T651 

1 

1 

0 

2 

2.5 

5 

l.E-07 

2 

15000 

YES 

CYCLE  Multiple 

18 

o 

11 

1 

0 

0 

1 

D7075-T651 

1 

1 

0 

2 

2.5 

5 

l.E-07 

2 

15000 

YES 

CYCLE  Multiple 

19 

o 

12 

1 

0 

0 

1 

D7075-T651 

1 

1 

0 

2 

2.5 

5 

l.E-07 

2 

15000 

YES 

CYCLE  Multiple 

20 

1 

12 

1 

0 

0 

1 

D7075-T651 

1 

1 

20000 

2 

2.5 

5 

l.E-07 

2 

15000 

YES 

CYCLE  Multiple 

21 

1 

12 

1 

0 

1 

1 

D7075-T651 

1 

1 

0 

2 

2.5 

5 

l.E-07 

2 

15000 

YES 

CYCLE  Multiple 

22 

1 

12 

1 

0 

0 

15 

D7075-T651 

1 

1 

0 

2 

2.5 

5 

l.E-07 

2 

15000 

YES 

CYCLE  Multiple 

23 

1 

12 

1 

0 

0 

1 

D7075-T651 

1 

12 

0 

2 

2.5 

5 

l.E-07 

2 

15000 

YES 

CYCLE  Multiple 

24 

1 

12 

1 

0 

0 

1 

D7075-T651 

1 

1 

0 

2 

2.5 

5 

l.E-07 

3 

054 

15000 

YES 

CYCLE  Multiple 

25 

1 

12 

1 

0 

0 

1 

D7075-T651 

1 

1 

0 

2 

2.5 

5 

l.E-07 

4 

05 

06 

15000 

YES 

CYCLE  Multiple 

26 

1 

12 

1 

0 

0 

1 

D7075-T651 

1 

1 

0 

2 

2.5 

5 

l.E-07 

5 

041 

15000 

YES 

CYCLE  Multiple 

27 

1 

12 

1 

0 

0 

1 

D7075-T651 

1 

1 

0 

2 

2.5 

5 

l.E-07 

6 

05 

15000 

YES 

CYCLE  Multiple 

28 

1 

12 

1 

0 

0 

1 

D7075-T651 

1 

1 

0 

2 

2.5 

5 

l.E-07 

7 

0432 

12900  15000 

YES 

CYCLE  Multiple 
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TableSO:  FAMSH  Test  Case  Fatigue  Life  Estimates 


Case 

Kn 

CGAP1.8  (FAMSH) 

FAMSH  1.50LINUX 

% 

difference 

1 

2.5 

98082.37 

98082.37 

0.00% 

5.0 

6341.44 

6341.44 

0.00% 

2 

2.5 

122287885.33 

122287885.33 

0.00% 

5.0 

4696715.79 

4696715.79 

0.00% 

3 

2.5 

28781.00 

28781.00 

0.00% 

5.0 

1860.58 

1860.58 

0.00% 

1 

2145999.00 

2145999.00 

0.00% 

10 

Notch  Strain  Exceeds  Fracture  Strain  For  Material 

4 

2.5 

233007031.05 

233007031.05 

0.00% 

5.0 

7181150.94 

7181150.94 

0.00% 

5 

2.5 

115019049.25 

115019049.25 

0.00% 

5.0 

3946915.39 

3946915.39 

0.00% 

6 

2.5 

221042834.00 

221042834.00 

0.00% 

5.0 

8424547.69 

8424547.69 

0.00% 

7 

2.5 

221042834.00 

221042834.00 

0.00% 

5.0 

8424547.69 

8424547.69 

0.00% 

8 

2.5 

221042834.00 

221042834.00 

0.00% 

5.0 

8424547.69 

8424547.69 

0.00% 

9 

2.5 

233007031.05 

233007031.05 

0.00% 

5.0 

7181150.94 

7181150.94 

0.00% 

10 

2.5 

115019049.25 

115019049.25 

0.00% 

5.0 

3946915.39 

3946915.39 

0.00% 

11 

2.5 

249151.21 

249151.21 

0.00% 

5.0 

10864.74 

10864.74 

0.00% 

12 

2.5 

220312.25 

220312.25 

0.00% 

5.0 

3457.82 

3457.82 

0.00% 

B 

2.5 

6.54 

6.54 

0.00% 

5.0 

0.42 

0.42 

0.00% 

14 

2.5 

122287885.33 

122287885.33 

0.00% 

5.0 

4696715.79 

4696715.79 

0.00% 

15 

2.5 

122287885.33 

122287885.33 

0.00% 

5.0 

4696715.79 

4696715.79 

0.00% 
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Case 

Kn 

CGAP1.8(FAMSH) 

FAMSH  1.50LINUX 

% 

difference 

16 

2.5 

122287885.33 

122287885.33 

0.00% 

5.0 

4696715.79 

4696715.79 

0.00% 

17 

2.5 

122287885.33 

122287885.33 

0.00% 

4696715.79 

4696715.79 

0.00% 

5.0 

18 

2.5 

122287885.33 

122287885.33 

0.00% 

5.0 

4696715.79 

4696715.79 

0.00% 

19 

2.5 

122287885.33 

122287885.33 

0.00% 

5.0 

4696715.79 

4696715.79 

0.00% 

20 

2.5 

48337584.70 

48337584.70 

0.00% 

5.0 

4563759.61 

4563759.61 

0.00% 

21 

2.5 

122287885.33 

122287885.33 

0.00% 

5.0 

4696715.79 

4696715.79 

0.00% 

22 

2.5 

40333647.35 

40333647.35 

0.00% 

5.0 

2148989.19 

2148989.19 

0.00% 

23 

2.5 

316442565.90 

316442565.90 

0.00% 

5.0 

10173569.74 

10173569.74 

0.00% 

24 

2.5 

28051631.66 

28051631.66 

0.00% 

5.0 

2624609.12 

2624609.12 

0.00% 

25 

2.5 

22761393.09 

2276B93.09 

0.00% 

5.0 

2007716.34 

2007716.34 

0.00% 

26 

2.5 

32232623256.35 

32232623256.35 

0.00% 

5.0 

968792708.66 

968792708.66 

0.00% 

27 

2.5 

32672040.32 

32672040.32 

0.00% 

5.0 

3145626.73 

3145626.73 

0.00% 

28 

2.5 

16197755.53 

16197755.53 

0.00% 

5.0 

2563178.05 

2563178.05 

0.00% 
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